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1  Overview  of  Mathematical  Accomplishments 

We  present  a  study  of  the  harmonic  analysis  of  finite  abelian  groups  A  and  finite 
groups  G  of  the  form  B  ,  the  semidirect  product  of  a  normal  abelian  group 
A  with  an  abelian  group  B.  The  theory  easily  generalizes  to  finite  groups  A<$  H} 
where  A  is  a  normal  abelian  group  and  if  is  an  arbitrary  group  whose  harmonic 
analysis  is  known.  In  particular,  it  includes  finite  groups  j4<$(E<$C),  where 
j4,  B  and  C  are  abelian  subgroups.  This  rich  class  of  finite  groups  contains  all 
crystallographic  groups. 

For  the  most  part  group  algebras  will  be  taken  over  the  complex  field  (split 
case),  but  a  section  on  group  algebras  of  abelian  groups  over  finite  fields  has  been 
included  as  an  indication  of  the  methods  that  must  be  developed  for  applications 
to  coding  theory. 

The  restriction  to  groups  G  =  A<§  B  and  their  generalizations  permit  har¬ 
monic  analysis  to  proceed  by  abelian  group  character  theory.  This  results  in 

•  a  direct  link  with  abelian  group  harmonic  analysis  and  consequently  with 
physical  interpretation. 

•  automatic  procedures  for  constructing  G-invariant,  irreducible  G-invariant 
and  direct  sum  decompositions  of  CG,  the  group  algebra  of  G  over  C,  into 
irreducible  G-invariant  subspaces. 

•  algorithms  for  computing  bases  of  irreducible  G-invariant  subspaces  and 
bases  compatible  with  direct  sum  decompositions  of  CG  into  G-invariant 
and  irreducible  G-invariant  subspaces  (spectral  bases). 

•  fastT  algorithms  for  representing  data,  i.e.,  elements  in  CG,  over  spectral 
basis. 

•  an  extensive  new  class  of  fast  unitary  transforms  for  data  analysis. 
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2  Introduction 


From  a  mathematical  perspective,  a  significant  part  of  classical  digital  signal 
processing  (DSP)  can  be  viewed  as  topics  in  finite  abelian  group  harmonic  anal¬ 
ysis  [  1,  6,  16,  4].  Fundamental  DSP  operations  such  as  convolution  and  the 
Fourier  transform  can  be  identified  with  group  algebra  multiplication  and  group 
algebra  direct  sum  decompositions  into  irreducible  group  invariant  subspaces. 
This  interplay,  often  implicit,  has  been  responsible  for  fast  algorithms  such  as 
the  FFT  [  7],  the  use  of  FFT  in  computing  large  size  convolutions  and  correla¬ 
tions,  and  more  recently  for  the  development  of  polynomial  transforms  [  27,  25] 
for  computing  convolutions. 

However  pleasing,  at  least  to  the  mathematically  inclined,  and  useful  this 
group  theoretic  approach  to  DSP,  it  is  based  on  a  seemingly  magical  relationship 
between  DSP  applications  and  finite  abelian  group  harmonic  analysis  and  the 
simplicity  of  this  harmonic  analysis.  This  magic  can  be  explained  by  the  power¬ 
ful  role  played  by  abelian  group  characters  in  providing  all  that  is  required  for 
finite  abelian  group  harmonic  analysis  and  the  physical  interpretation  of  these 
characters  as  frequency. 

For  some  mathematicians,  this  group  interpretation  of  DSP  has  raised  the 
potential  of  an  equally  important  application  of  nonabelian  group  harmonic 
analysis  to  DSP,  especially  in  the  construction  of  group  transforms  and  group 
filers  generalizing  the  classical  Fourier  transform  and  convolutional  filters.  The 
works  of  R.  Holmes  [  14,  15],  M.  Karpovski  and  E.  Trachtenberg  [  18]  are  the 
basis  of  much  of  the  research  in  this  direction.  Similar  ideas  in  coding  theory 
have  been  introduced  by  F.J.  MacWilliams  [  21].  These  efforts  have  shown 
some  promise  but  for  the  most  part  are  more  interesting  to  mathematicians 
than  DSP  engineers.  A  more  successful  application  has  been  to  fast  algorithm 
design  [  22,  23,  2,  3,  10]. 

During  the  last  ten  years  considerable  effort  has  taken  place  to  extend  the 
success  and  range  of  applicability  of  nonabelian  group  methods  to  the  design  of 
new  filters  and  spectral  analysis  methodologies  [  8,  9,  5],  as  an  image  processing 
tool  [  11,  13,  19,  20]  and  most  recently  as  an  image  processing  tool  combined 
with  graph  theoretic  modeling  of  image  data  [  12,  24]. 

Many  efforts  at  finding  a  significant  role  for  nonabelian  group  theory  in  DSP 
and  imaging  applications  as  well  as  in  coding  and  communication  theory  have 
been  limited  by  some  or  all  of  the  following. 

•  The  choice  of  an  appropriate  group  or  groups  in  a  given  application  is  not 
obvious.  Often  the  groups  considered  are  those  which  are  best  known  to 
the  researcher. 

•  The  lack  of  a  large  class  of  groups  whose  harmonic  analysis  is  sufficiently 
understood  for  meaningful  applications.  Often  the  dihedral  group  with  or 
without  justification  is  the  test  example. 
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•  The  need  to  develop  a  conceptual  framework  which  relates  group  harmonic 
analysis  to  physically  interpretable  results. 

•  Fast  algorithms  may  exist  but  are  difficult  to  code  as  they  do  not  easily 
lend  themselves  to  modification. 

•  Relevant  models  for  applications  which  are  not  only  a  rephrasing  of  known 
methods  or  what  is  immediately  available. 

In  the  following  three  sections,  we  present  a  detailed  account  of  finite  abelian 
group  harmonic  analysis.  This  theory  is  well-known  to  most  mathematicians  but 
may  not  be  familiar  to  DSP  engineers.  It  is  the  starting  point  for  the  nonabelian 
group  harmonic  analysis  developed  in  this  work  since  as  we  will  see  both  the 
abelian  group  and  nonabelian  group  theories  share  many  common  features. 

The  finite  nonabelian  groups  in  this  study  have  the  form  G  =  A<$B,  the 
semidirect  product  of  a  normal  abelian  group  A  with  an  abelian  group  B .  Re¬ 
sults  easily  extend  to  finite  groups  where  H  is  an  arbitrary  group  whose 

harmonic  analysis  is  known.  This  includes  finite  groups  of  the  form  A<$  (B<$  G), 
where  A,  B  and  C  are  abelian  groups  and  in  particular  all  crystallographic 
groups.  In  the  applications  part,  we  will  see  how  these  nonabelian  groups  pro¬ 
vide  for  new  classes  of  imaging  and  DSP  models. 

The  language  of  group  algebras  usually  over  the  complex  field  will  be  used 
throughout.  Some  readers  may  be  more  familiar  with  the  identification  of  the 
group  algebra  C G  with  the  space  of  all  complex  valued  functions  on  G  under 
G-convolution.  Harmonic  analysis  over  G  usually  includes  a  description  of  the 
(left)  G-invariant  and  irreducible  G-invariant  subspaces  and  direct  sum  decom¬ 
positions^  CG  into  irreducible  G-invariant  subspaces. 

For  an  abelian  group  A,  complete  answers  to  these  problems  can  be  given 
in  terms  of  the  characters  of  A.  Moreover,  these  descriptions  are  especially 
simple.  For  example  the  irreducible  A-invariant  subspaces  coincide  with  the 
one-dimensional  subspaces  spanned  by  the  characters.  The  characters  of  A 
determine  an  orthogonal  basis  of  CA  and  the  finite  Fourier  transform  relates  the 
delta  basis  A  of  C  A  with  the  basis  of  characters  of  A.  This  is  an  extremely  nice 
answer  for  applications  since  the  characters  of  A  can  be  physically  interpreted 
as  frequencies. 

For  finite  groups  A<$  B,  the  characters  of  A  and  B  play  an  equally  important 
but  more  complicated  role  in  the  harmonic  analysis  of  A<$B.  This  theory  is 
developed  in  section  6  and  extensions  to  A<$  (j3<$  C)  and  beyond  will  be  devel¬ 
oped  in  time.  We  will  have  as  a  result  a  large  class  of  nonabelian  groups  with 
the  potential  of  wide  applicability  to  DSP  and  imaging  whose  harmonic  analysis 
is  known  in  detail.  This  harmonic  analysis  is  based  on  abelian  group  characters 
which  directly  link  this  harmonic  analysis  to  familiar  DSP  and  imaging  concepts. 
Moreover  we  have  developed 

•  algorithms  for  computing  bases  of  irreducible  G-invariant  subspaces  and 
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bases  compatible  with  direct  sum  decompositions  of  C G  into  G-invariant 
and  irreducible  G-invariant  subspaces  (spectral  bases). 

•  fast  algorithms  for  representing  data,  i.e.,  elements  in  CG  over  spectral 
bases. 

•  an  extensive  new  class  of  fast  unitary  transforms  for  data  analysis. 

The  direct  sum  decompositions  of  CG  studied  in  this  work  are  closely  related 
to  but  are  significantly  different  from  those  leading  up  to  the  Fourier  transform 
of  CG  [  10,  26].  Our  decompositions  are  finer  but  are  not  generally  uniquely 
determined. 

Throughout  the  mathematics  part  A  and  B  denote  finite  abelian  groups  with 
composition  given  by  multiplication.  The  identity  is  always  denoted  by  1  and 
the  inverse  of  x  €  A  by  x-1.  The  order  of  a  set  is  the  number  of  elements  in 
the  set. 

For  much  of  the  general  theory  K  denotes  a  field,  finite  or  infinite,  whose 
characteristic  does  not  divide  the  orders  of  A  and  B.  However  for  DSP  and 
imaging  applications,  K  is  the  field  of  complex  numbers.  A  section  on  group 
algebras  of  A  over  finite  fields  is  included  as  an  indication  of  the  methods  used 
in  applications  to  coding  theory. 
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3  Characters  of  Finite  Abelian  Groups 

A  mapping  r  :  A  — >  Kx  is  called  a  character  of  A  over  K  if  r  is  a  homomorphism 
of  the  group  A  into  the  multiplicative  group  Kx  of  nonzero  elements  in  K , 

r(xy)  =  r(x)r(y),  x,  y  €  A. 

Denote  by  ch(A,  A")  the  set  of  all  characters  of  A  over  K. 

For  each  a:  €  A,  denote  by  gp^(x)  the  group  generated  by  x  in  A, 

gPyl(*)  =  {!»*,  •••»  z*-1}, 

where  the  order  R  of  gpA(x)  is  the  smallest  positive  integer  satisfying  xR  =  1. 
We  call  R  the  order  of  x. 

If  r  is  a  character  of  A  over  K  and  x  €  A,  then  r(x)  completely  determines 
the  values  of  r  on  #pa(x)  by 

r(x!)  =  t(x)*,  /gZ. 

In  particular,  xL  —  1  implies  r(x)  is  an  L-th  root  of  unity  in  K . 

Denote  by  Un(K)  the  multiplicative  group  of  all  N-th  roots  of  unity  in 
K.  The  order  of  Un{K)  divides  AT.  Since  every  x  G  A  satisfies  xN  =  1, 
r(x)  e  Uw(K),  for  every  character  r  of  A  over  K . 

A  spfc'ts  over  A"  if  the  number  of  characters  of  A  over  K  is  the  order  of  A. 
We  will  eventually  show  that  this  is  the  maximal  number  of  possible  characters 
of  A  over  K . 

In  the  following  examples  Cjv  is  the  cyclic  group  of  order  N  having  generator 
x. 

Example  1  There  exists  N  characters  of  Cw  over  C  defined  by 
rn(x)  =  e27rtft,  0  <  n  <  N. 

The  multiplicative  group  Un(K)  is  a  cyclic  subgroup  of  Kx  whose  order  L 
divides  N.  If  7  €  then  7L  =  1.  The  smallest  positive  power  of  7  equal 

to  1  is  called  the  order  of  7.  Denoting  a  generator  of  Un(K)  by  a,  we  have  that 

UN(K)  =  {l,a,...,aL~1} 

and  a  has  order  L.  We  say  that  a  is  a  primitive  N- th  root  of  unity  if  the  order 
of  a  is  N.  K  has  a  primitive  N- th  root  of  unity  if  and  only  if  there  exists  exactly 
N  N- th  roots  of  unity  in  K. 

Cn  splits  over  K  if  and  only  if  K  has  a  primitive  N~ th  root  of  unity  a.  The 
N  characters  of  over  K  are  defined  by 

rn(x)  =  an,  0  <  n  <  N. 
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Example  2  cyclotomics 

Suppose  K  =  GF(pR ),  the  field  of  order  pR,  p  a  prime.  Since  Kx  is  a 
cyclic  multiplicative  group  of  order  pR  —  1,  a  generator  a  of  Kx  is  a  primitive 
(pR  —  l)-th  root  of  unity  and  the  order  of  any  power  of  a  must  divide  pR  -  1. 
If  N  divides  pR  —  1,  then 

7  =  aM,  pR  -  1  =  MN, 

is  a  primitive  N- th  root  of  unity  in  Kx.  We  see  that  CN  splits  over  GF(pR)  if 
and  only  if 

pR  =  1  mod  N. 

In  particular,  if  C n  splits  over  GF(pR)  then  p  does  not  divide  N. 

Example  3  The  finite  field  K  =  GF(27)  has  a  primitive  26-th  root  of  unity  a, 

K*  ={1,  a,  a25} 

with  a26  the  smallest  positive  power  of  a  equal  to  1. 

Example  4  Continuing  the  preceding  example,  a 2  has  order  13,  a13  has  order 
2  and  ar  has  order  26  whenever  r  is  relatively  prime  to  26. 

Example  5  Ce  and  C%  split  over  GF{ 25). 

Example  6  C4  splits  over  GF(  9). 

In  general,  the  number  of  characters  of  Cn  over  K  is  equal  to  the  number  of 
iV-th  roots  of  unity  in  K.  If  Un(K)  has  order  L  and  a  is  a  generator  of  Un(K), 
the  L  characters  of  Cn  over  K  are  defined  by 

ri(x)  =  qz,  0  <1  <  L. 

If  K  does  not  contain  a  primitive  N-th  root  of  unity,  then  the  number  of  char¬ 
acters  L  of  Cn  over  K  is  less  than  N.  In  fact,  L  is  a  proper  divisor  of  N. 

Example  7  R  has  two  roots  of  unity,  1  and  —1.  If  N  is  odd,  1  is  the  only  7V-th 
root  of  unity  in  R  and  C n  has  exactly  one  character  over  R.  If  N  is  even,  1  and 
—  1  are  the  only  AT-th  roots  of  unity  in  R  and  Cn  has  exactly  two  characters 
over  R. 

Example  8  GF( 5)  has  exactly  two  2-th  roots  of  unity,  1  and  4,  and  four  4-th 
roots  of  unity,  1,  2,  3,  4.  There  are  exactly  two  6-th  roots  of  unity  1  and  4, 
and  exactly  two  characters  of  Cq  over  GF( 5).  There  are  exactly  four  8-th  roots 
of  unity,  1,  2,  3,  4  and  exactly  four  characters  of  Cg  over  GF( 5).  In  general  if 
N  =  2  mod  4,  there  are  exactly  two  characters  of  Cn  over  GF(5)  and  if  N  =  0 
mod  4,  there  are  exactly  four  characters  of  Cn  over  GF( 5). 
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For  K  =  Q  or  K  =  GF{pR))  p  a  prime  not  dividing  JV,  the  condition  that 
Cn  splits  over  K  places  severe  conditions  on  N.  One  solution  is  to  construct 
a  minimal  field  extension  E  of  K  over  which  Cn  splits.  E  contains  K  as  a 
subfield,  has  a  primitive  N- th  root  of  unity  and  there  exists  no  proper  subfield 
of  E  containing  K  having  a  primitive  iV-th  root  of  unity.  There  is  an  elaborate 
theory  for  construction  such  a  field  extension  over  an  arbitrary  field  K  whose 
characteristic  p  does  not  divide  N.  The  resulting  field  extension  E  is  called  the 
splitting  field  of  the  polynomial  xN  —  1  over  K  and  is  an  example  of  a  Galois 
extension.  We  will  say  more  about  Galois  extensions  over  Q  and  GF(pR)  in 
chapter  4. 

For  the  fields  Q  and  GF(pR),  the  construction  of  minimal  field  extensions 
over  which  Cn  splits  is  simple  and  we  will  restrict  our  attention  to  these  cases. 

For  Q,  the  cyclotomic  field  Q (£#),  =  e2™^,  is  the  minimal  field  exten¬ 

sion  of  Q  over  which  Cn  splits.  For  GF(pR)  the  construction  is  slightly  more 
involved.  For  simplicity  we  consider  the  case  GF(p ),  where  p  does  not  divide 
N. 

Viewing  p  as  an  element  in  Z/AT,  since  p  and  N  are  relatively  prime,  p  is 
in  the  multiplicative  group  U(N)  of  units  of  Z/N  and  some  power  of  p  equals 
1  mod  N.  The  smallest  such  power  it!,  called  the  order  of  p  in  U(N ),  defines 
the  minimal  extension  field  GF(pR)  over  which  Cn  splits.  In  general,  for  any 
integer  s  >  0,  there  exists  a  minimal  extensional  field  of  GF(ps)  over  which  Cat 
splits. 

Example  9  GF( 25)  is  the  minimal  extension  field  of  GF( 5)  over  which  Ce  and 
C&  splits.  GF{ 9)  is  the  minimal  extensional  field  of  C?F(3)  over  which  C4  splits. 

Consider  the  direct  product  Cnx  x  Cn2  of  cyclic  groups  Cnx  and  Cn2-  De¬ 
noting  generators  of  Cnx  and  Cn 2  by  X\  and  x2,  every  element  in  C^1  x  C^2 
can  be  written  uniquely  in  the  form 

xTx22'  0<m<iV1>0<n2<  N2. 

The  characters  of  Cnx  x  Cn2  over  K  ar€  completely  determined  by  the  characters 
of  Cni  and  Cn2  over  K*  If  t  is  a  character  of  Cnx  x  Cn2  over  we  can  define 
characters  T\  and  r2  of  C^i  and  Cn2  over  K  by 

ri(rEj)  =  r(x),  r2(x2)  =  r(x2). 

r  is  completely  determined  by  r\  and  r2  by 

r(xilx22)  =  Ti(xx)nir2(x2)n2i  0  <  ni  <  Ni,  0  <  n2  <  7V2, 

and  we  can  write  r  =  t\  0  r2.  Conversely,  if  T\  and  r2  are  characters  of  Cnx 
and  Cn2  over  K,  then  r\  <8>  r2  is  a  character  of  C^x  x  Cn2  over  K. 

Cni  x  Cn2  splits  over  K  if  and  only  if  Cnx  and  Cn2  split  over  K,  This 
will  be  the  case  if  and  only  if  K  has  a  primitive  N\- th  root  of  unity  ax  and  a 
primitive  N2- th  root  of  unity  ot2 . 
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Example  10  The  iVi  JV2  characters 

rn,  n  =  (ni,n2),  0  <  n\  <  N\,  0  <  n2  <  N2 , 
of  Cn1  xCjv2  over  C  are  defined  by 

Tn(x^x^)  =  eiwiS^e9wiS^‘,  0  <  mj  <  Nu  0  <  m2  <  JVa. 

If  If  has  primitive  iN^-th  and  AT2-th  roots  of  unity  oti  and  a2,  then  the  N1N2 
characters 

rn,  n  =  (ni,n2),  0  <  rii  <  Ni,  0  <  n2  <  AT2, 
of  x  Cn2  over  if  are  defined  by 

TmipT**2)  =  a?imio£2Tn2,  0  <  mi  <  ATi,  0  <  m2  <  AT2. 


Denote  the  least  common  multiple  of  N\  and  N2  by  [iVx ,  iV2 ] .  a\a2  is  a 
primitive  [JVi,  AT2]-th  root  of  unity.  Generally  if  K  has  a  primitive  [JVi,iV2]-th 
root  of  unity  a,  then 

ai  =  aLl,  a2  =  aLs,  [iVi , iVa]  =  L1JV1  =  L2N2 

are  primitive  Ni- th  and  AT2-th  roots  of  unity  in  K. 

K  has  primitive  A^-th  and  iV2-th  roots  of  unity  if  and  only  if  K  has  a 
primitive  [iVi,  AT2]-th  root  of  unity.  Cn1  x  Cn2  splits  over  K  if  and  only  if  K 
has  a  primitive  [Ni ,  iV2]-th  root  of  unity. 

We -can  state  the  condition  for  Cn1  x  Gn2  to  split  over  K  completely  in  terms 
of  the  abelian  group  defined  by  Cnx  x  Cpi2  and  not  on  its  representation.  Since 
the  order  of  each  element  in  Cnx  divides  N\  and  the  order  of  each  element  in 
Cn2  divides  AT2,  [iVi,  N2]  is  the  least  common  multiple  over  all  orders  of  elements 
in  CNl  x  CN2. 

In  general  the  number  of  characters  of  Cnx  xCn2  over  K  is  the  product  of 
the  number  of  characters  of  Cnx  over  K  with  the  number  of  characters  of  Cjv2 
over  K.  This  number  is  equal  to  L\L2  where  L\  is  the  number  of  Ni-th  roots 
of  unity  and  L2  is  the  number  of  AT2-th  roots  of  unity  in  K.  If  a\  and  a2  are 
generators  of  Unx(K)  and  Upj2{K))  the  L\L2  characters  of  Cnx  x  Cn2  over  K 

Th  1  =  (^1  >  ^2)5  0  <  l\  <  Li,  0  <  l2  <  L2, 


are  given  by 

lx?2)  =  a[imialr?m\  0  <  mi  <  Nu  0  <  m2  <  N2. 

The  results  of  this  section  easily  extend  to  an  arbitrary  finite  number  of 
cyclic  groups 

A  —  Cpi j  x  •  ■  •  x  Cnr. 
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Each  character  r  of  A  over  if  is  uniquely  represented  by 

r  =  ti  ®  •  •  *  0  tj*, 

where  rr  is  a  character  of  CVr  over  if ,  1  <  r  <  i2.  A  splits  over  if  if  and  only 
if  Cjvr  splits  over  if ,  1  <  r  <  R  and  in  general  the  number  of  characters  of  A 
over  if  is  equal  to  L  =  Li  •  •  ■  L;??  where  Lr  is  the  number  of  characters  of  Cn r 
over  if ,  1  <  r  <  R.  Since  every  finite  abelian  group  A  is  the  product  of  a  finite 
number  of  finite  cyclic  groups,  this  completes  the  description  of  Ch(A  :  K)  for 
every  finite  abelian  group  A  and  field  if. 

4  Group  algebra  of  A 

The  group  algebra  K A  of  A  over  K  is  the  if -vector  space  of  all  formal  sums 

/  =  £  /(*)*,  /(*)  €  K, 

x€A 


under  . 

/  + 9  =  X)  (/(*)  +  S(x))  *>  f,g€KA, 

xsa 

«/= E  (a/(®))x,  a  Gif,  /Gif  A, 

with  if -algebra  multiplication 

/ff  =  (  E  fix)9{x~ly) )  2/.  /«  5  €  if  A 

\x€A  / 

The  additions  and  multiplications  inside  the  brackets  are  taken  in  if.  Since  A 
is  abelian,  the  if-algebra  multiplication  is  commutative. 

Denote  by  L(A\K),  the  if-vect or  space  of  all  if-valued  functions  on  A. 
Every  /  G  L(A  :  if)  defines  a  formal  sum  in  if  A 

f  =  E  /(*)* 

and  we  can  identify  the  if- vector  space  if  A  with  the  if- vector  space  L(A\  if ). 
The  multiplication  /  •  g  in  if  A  corresponds  to  the  standard  convolution  in 
L(A;K) 


f  *  9(y)  =  ^2  f(x)g(x  1y),  y  e  A,  f,  g  e  L{A;K). 


Under  this  identification  the  delta  function  6y,  y  £  A  corresponds  to  the  formal 
sum  in  KA  __ 

Y  Mx)x 

x€A 

which  we  denote  by  y.  In  this  way,  we  can  view  A  as  a  subset  of  KA.  A  is  a 
basis  of  the  Af-vector  space  KA  and  KA  has  dimension  AT,  the  order  of  A.  The 
identity  1  in  A  is  the  identity  of  the  AT-algebra  KA. 

We  will  usually  use  the  term  basis  for  any  subset  of  vector  space  which  if 
ordered  is  a  basis  in  the  usual  sense.  In  summation  expressions,  there  is  no 
loss  in  doing  so.  However  whenever  matrices  are  involved,  an  ordering  must  be 
specified  even  if  implicitly. 

For  y  £  A  and  /  £  KA 

vf  =  Y  f(x)yx  =  Y  f(y~lx)x- 

x£A 

We  call  yf  the  translation  of  /  by  y.  The  term  f(y~lx)  in  the  right-hand 
summation  identifies  yf  with  the  usual  definition  of  translation  in  engineering 
terminology.  For  y  £  A,  the  operator  L(y)  of  KA  defined  by 

L(y)f  =  yf,  f£  KA, 

is  a  linear  isomorphism  of  the  if-vector  space  KA  and  we  have 
L(zy)  =  L{z)L{y ),  z,y  £  A 

L(z~x)  =  L(z)~l ,  z  €  A 

where  L(z)L(y)  and  Li^z)"1  denote  composition  and  inversion  of  operators  in 
KA.  Observe  that  for  every  y  £  A,  L(y)  acts  as  a  permutation  on  the  basis  A 
of  KA. 

Example  11  Ordering  C4  by  successive  powers  of  a  generator  x, 

1,  x,x2,x3, 

the  matrix  of  L(xm)  relative  to  the  resulting  basis  C4  of  KC4  is  S™,  0  <  m  <  4, 
where  54  is  the  4-point  cyclic  shift  matrix. 

In  general,  ordering  Cn  by  successive  powers  of  a  generator  x ,  the  matrix  of 
L(xm)  relative  to  the  resulting  basis  C x  of  KCn  is  S^} ,  0  <  m  <  N,  where  Sn 
is  the  AT-point  cyclic  shift  matrix. 

Example  12  Ordering  C4  XC3  by  first  ordering  C4  and  C3  by  successive  powers 
of  generators  X\  and  X2  and  then  by  imposing  the  lexicographic  ordering  on 

C4  x  C3 

1,  £2,  %\\  X\X2,  X\x\\  x\ ,  x\x2 ,  x\x\\  x\ ,  x\x2 ,  x\x\, 

the  matrix  of  Lfe^'x™2)  relative  to  the  resulting  basis  C4  x  C3  of  K(C4  x  C3) 
is  S™1  0  S™2 ,  0  <  mi  <  4,  0  <  m2  <  3. 
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In  general  ordering  the  direct  product 


CNlx  -x  CNr 

by  first  ordering  the  factors  C^r  by  successive  powers  of  generators  xrj  1  <  r  < 
R)  and  then  by  imposing  the  lexicographic  ordering  on  the  direct  product,  the 
matrix  of 

•••*;*) 

relative  to  the  resulting  basis  Cnx  x  •  ■  •  x  C^R  of  K(Cnx  x  •  •  •  x  C^r)  is 

Denote  by  GL{KA\  if)  the  group  of  linear  isomorphisms  of  the  if-ve ctor 
space  if  A.  The  mapping 

L:A->GL{KA]K) 

is  an  isomorphism  of  the  group  A  into  the  group  GL(KA,  if).  We  say  that  L 
is  a  representation  of  A  on  KA  and  call  L  the  regular  representation  of  A. 

A  subspace  V  of  the  if -vector  space  KA  is  called  A-invariant  if  for  each 
ye  A 

yV  =  {yf  :  /  €  V}  C  V. 

If  V  is  A-invariant  and  /(x),  x  €  A,  is  a  coefficient  set  for  a  formal  sum  in  V) 
then  for  each  y  £  A,  f(y~xx)^  x  €  A,  is  also  a  coefficient  set  for  a  formal  sum 
in  V.  A-invariant  subspaces  of  KA  can  be  identified  with  translation-invariant 
subspaces  of  if-valued  functions  of  A. 

Example  13  The  subspace  of  if  A  spanned  by  the  formal  sum 

x£A 

is  the  A-invariant  subspace  of  if  A  consisting  of  all  formal  sums  in  if  A  unchanged 
by  the  action  of  the  operators  L(y)}  y  €  A. 

For  a  subgroup  B  of  A,  we  say  that  a  subspace  V  of  the  if-vector  space  if  A 
is  B-invariant  if  for  all  y  €  B, 

yV  C  V. 

The  ^-invariant  subspace  of  if  A 

{/  e  KA  :  yf  =  /,  for  all  y  e  B} 

can  be  identified  with  the  subspace  of  B-periodic  if -valued  functions  on  A. 

For  g  €  if  A  define  the  operator  L(g)  of  if  A  by 

L(g)f  =  gf,  f  e  KA. 
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Since 


L{g)f  =  gf  =  J2  9(y)(yf )  =  Y  9{v)Hv)f  = 

y€A  y€A 

we  have  that  _ 

L(y)  =  Y  9{y)L(y) 

y€A 

linearly  extends  the  domain  of  definition  of  L  from  A  to  KA.  In  particular,  the 
two  definitions  coincide  on  A . 

L(g) ,  g  G  KA)  is  a  homomorphism  of  the  if -vector  space  KA  but  is  not 
necessarily  an  isomorphism  as  is  the  case  with  L(y)>  y  G  A.  We  can  have 
gf  =  0  for  /,  g  G  KA  with  both  /  and  g  not  zero. 

Example  14  Continuing  example  11,  the  matrix  of  L(g)  g  G  K(C4 ),  relative 
to  the  basis  C4  is 

Ci{g)=YJg{xm)S?. 

m— 0 

Defining  g  G  KA  by 

g  =  b(xm)]0<m<4 

we  have 

C4(g)  =  [  g  S4g  S4g  S43g]. 

If  we  take  g  =  2/»  ^hen  C4G7)  1S  the  singular  4x4  matrix  of  all  ones. 

In  general,  ordering  Cn  by  successive  powers  of  a  generator  a:,  the  matrix  of 
L(g))  g  G  KCn,  relative  to  the  resulting  basis  Cn  is 

CN(g)=Y,g(xm)S%. 

m=0 

Defining  g  G  KN  by 

g  =  [s(*m)]0<m<AM 

we  have 

cJV(s)  =  [g  sN g  •••  S%-' g]. 

Linear  combinations  over  if  of  powers  of  Sjv  are  called  iV-point  circulant 
matrices  over  K.  Such  matrices  are  completely  determined  by  their  0-th  column. 

Cat  can  be  realized  relative  to  the  basis  Cw  as  the  matrix  product  of  the 
circulant  matrix  Cw(g)  with  the  vector  f  formed  by  ordering  the  values  of  / 
according  to  successive  powers  of  the  generator. 


Y  y(y)L{y)  f 

ySA  ) 
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Example  15  Continuing  example  14,  for  /  and  g  in  KC 4,  the  product  h  =  fg 
is  given  by 

h  =  C4(ff)f. 

Example  16  Continuing  example  12,  the  matrix  of  L(g ),  p  €  if(C*4  x  C3) 
relative  to  the  basis  C4  x  C3  is 

E  E 

mi=0  7712=0 

Denote  by  C(m)  the  3x3  circulant  matrix  having  0-th  column 

g(x?) 

g(x™x 2)  ,  0  <  m  <  4. 

.  . 

The  matrix  of  L(g)  can  be  written  as 

'  c(0)  c(3)  c(2)  c(l)  ' 

c(l)  c(0)  c(3)  c(  2) 

c(2)  c(l)  c(0)  c(3) 

.  c(3)  c(2)  c(l)  c(0)  . 

which  is  an  example  of  a  block- circulant  matrix  having  circulant  blocks . 

Denote  by  gl(KA]K)  the  K- algebra  of  all  homomorphisms  of  the  K -vector 
space  KA.  The  mapping 

L:KA-+gl{KA;K) 

is  a  if-algebra  isomorphism  of  KA  into  gl{KA\K). 

For  an  A-invariant  subspace  V  of  KA  and  /  €  V',  we  have  for  all  g  €  KA , 

L{g)f  =  E  9(y)Hy)f  e  v 

and  y  is  KA- invariant.  The  equivalence  of  A-invariance  and  if  A-invariance 
will  be  used  throughout  this  work. 

A  homomorphism  T  of  the  if -vector  space  KA  is  called  A-invariant  if  for 
all  y  €  A, 

L(y)T  =  TL(y). 

Example  17  L(g),  g  €  KA  is  A-invariant. 
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In  fact,  the  operators  L(g ),  g  €  KA  are  the  only  ^-invariant  homomorphisms 
of  KA.  For  an  ^-invariant  T  and  f  £  KA 

yT(f)  =  T(yf),  y  €  A. 

In  particular,  with  /  =  1  in  KA, 


T(y)  =  T(l)y 


and 


tu)  =  E  f(y)T(y)  =  T(!)  E  /(»)*  =  to/- 

The  claim  is  proved  since 

T  =  L(g),  g  =  T(  1). 


5  Fourier  transform  over  A 

Decompositions  of  KA  into  direct  sums  of  ^-invariant  subspaces  play  two  im¬ 
portant  roles:  algorithm  design  for  computing  products  in  KA  and  spectral 
analysis.  The  Fourier  transform  over  A  is  the  simplest  and  most  frequently 
occurring  example. 

A  character  r  of  A  over  K  determines  the  formal  sum  in  KA 

r  =  ^  r(x)x. 

xeA 

Multiplication  of  characters  will  always  be  taken  in  KA  with  the  warning  that 
in  many  places,  another  multiplication  is  defined  under  which  ch {A\K)  is  a 
group.  The  importance  of  characters  in  studying  A-invariant  subspaces  of  KA 
is  contained  in  the  following  result. 

Theorem  1  For  f  €  KA  and  r  a  character  of  A  over  K, 

fr  =  /(r)r, 

where  f(r)  €  K  is  given  by 

Hr)  =  E  /(vMir1)- 

y€A 


Proof  For  y  €  A 

yT  =  e  T(x)yx  =  E  t(v"1*)* = E  T^x 

x€A  xEA  xEA 

=  r{y~l)T. 
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The  theorem  follows  from 


fr  =  53  f(y)(yT) = 

yeA 


53  /(2/)r(2/  *)  T' 

yeA  j 


For  two  characters  r  and  A  of  A  over  if,  since  rA  =  At,  we  have 

rA  =  aA  =  ar  =  At, 

where  a  =  t(A)  =  A(t).  If  t  ^  A,  then  a  =  0  and  tA  =  0.  If  t  =  A,  then 

a  =  53  r(y)T(y-1)  =  5Z 1  = 

and  t2  —  Nt,  proving  the  following. 

Corollary  1  For  two  characters  r  and  X  of  A  over  if , 


tA 


_  f  JVr,  t 

"l  o,  T 


=  A, 
#A. 


In  the  following  discussion,  we  abbreviate  ch(A;  if )  by  A*  and  view  A*  as  a 
subset  of  if  A. 

Corollary  2  A*  is  a  linearly  independent  subset  in  the  if -vector  space  KA. 
Proof  Suppose 

0=  y,  Q(r)r’  a(r) G K- 

t€A* 

By  corollary  1,  for  any  A  €  A* 


X  ]T  o(t)t  =  Na( X)X  =  0 

r£A* 

and  a(A)  =  0.  Since  this  holds  for  any  A  €  A*,  the  corollary  follows. 

By  the  theorem,  the  if -subspace  spanned  by  a  character  t  over  if  is  A- 
invariant 

if  At  =  if  t. 


By  corollary  2,  we  have  the  direct  sum  of  one-dimensional  A-invariant  subspaces 
of  if  A, 

53®tfr. 

re  a* 
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5.1  Split  case 

Suppose  A  splits  over  K .  By  corollary  2,  A*  is  a  basis  of  the  K-vectov  space 
KA  and  we  have  the  following  result. 

Theorem  2  KA  is  the  direct  sum  of  one- dimensional  A-invariant  subspaces 

KA=  ®Kt • 

rGA * 


By  theorem  2, 

1  =  ^  a(r)r,  a(T)  e  K. 

t€A* 

For  any  A  G  A*, 

A  =  A  •  1  =  ^  a(r)Ar  —  7Va(A)A 

rem¬ 
and  a(A)  =  proving  the  following. 

Corollary  3 


-  T 

Nrh 


T. 


By  corollary  3, 

/  =  ^E-fr  =  ^E  /(r)T 

t£A*  reA* 

is  the  expansion  of  /  over  the  basis  A*  of  KA.  The  coefficient  set  of  this 
expansion  (up  to  scale  multiple 

/>),  TSA\ 

is  called  the  Fourier  transform  of  /  in  KA.  By  theorem  1, 

fir)  =  53 /(vMv"1)- 

y6  A 

The  linear  isomorphism  that  maps  the  coefficient  set  of  the  expansion  of  /  over 
the  basis  A 

fix),  xeA, 

onto  the  coefficient  set  of  the  expansion  of  /  over  the  basis  A*  (up  to  scale 
multiple  jr) 

fir),  t  £  A* 


18 


is  called  the  Fourier  transform  over  KA .  The  inversion  that  maps  coefficient 
sets  over  A*  onto  coefficient  sets  over  A 

/(*)  -  jj  5Z  /Vm*).  1  €  A, 

r€A * 

is  called  the  inverse  Fourier  transform  over  KA. 

Example  18  Examples 

The  basis  A*  of  KA  is  a  diagonalizing  basis  for  the  operators  L(g),  g  £  KA, 
since 

L(q)t  =  gr  =  g(r)r ,  r  €  A* 
and  a  diagonalizing  basis  for  multiplication  in  KA. 

Theorem  3  For  f  and  g  in  K A, 

f9  =  JjlL  /(r)»(T)r- 

reA* 

Proof 

/s  -  HsHfeH 

—  =  Jp  E  E-/»^(A)r a 

=  ^  E  f(T)9(r)T, 

reA* 

completing  the  proof. 

Corollary  4  ///  €  XA  satisfies  f2  =  0,  then  f  =  0. 

Proof  By  theorem  3 

/!  =  5S/|r|V  =  t 

implying,  since  A*  is  a  basis  of  iCA,  that 

/(r)  =  0,  re  A*, 
and  /  =  0  proving  the  corollary. 

The  theorem  leads  to  an  algorithm  for  computing  fg,  f,g€  KA. 
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•  Compute  the  Fourier  transform 

/». 

•  Compute  the  Fourier  transform 

9{t), 

•  Compute  the  product  in  K 

f{T)g(T),  re  A*. 

•  Compute  the  inverse  Fourier  transform  of  these  products 

/ff. 

The  importance  of  this  algorithm  to  compute  fg ,  /,  g  G  ATA,  based  on  the 
existence  of  fast  algorithms  to  compute  the  Fourier  transform. 

An  A-invariant  subspace  W  of  K A  is  called  irreducible  if  the  only  A-invariant 
subspaces  of  IF  are  (0)  and  W.  For  r  G  A*,  Kr  is  an  irreducible  A-invariant 
subspace  of  KA.  In  the  split  case  the  A- invariant  and  irreducible  A- invariant 
subspaces  of  KA  can  be  completely  described  in  terms  of  the  subspaces  Kt, 
r  €  A*. 

Theorem  4  If  W  is  an  A-invariant  subspace  of  KA,  then 

W=Y,  ®Kt  =  KAe, 

_ _  r€  A 

where  A  =  W  n  A*  and  e  =  j,  J2T eA  T- 
Proof  For  /  £W  and  A  €  A* 

»/4E  hr)Xr  =  i/(A)A. 

Since  IF  is  A-invariant, 

/( A)A  G  IF,  A  €  A*. 

If  /(A)  ^  0,  then  A  G  IF  proving  that  every  /  G  IF  is  contained  in  the  if-linear 
span  of  the  set  of  characters  contained  in  IF,  proving  IF  =  £r€A 

Since  e  €  W  and  W  is  A-invariant,  KAe  C  IF.  For  any  A  G  A,  Ae  =  ^A  € 
KAe ,  proving  IF  C  i^Ae,  completing  the  proof  of  the  theorem. 

The  factor  ^  in  the  definition  of  e  has  been  chosen  so  that  e2  =  e,  the 
relevance  of  which  will  be  made  clear  in  the  next  section. 

If  IF  is  an  irreducible  A-invariant  subspace  of  KA,  then  by  theorem  4,  IF 
contains  a  unique  r  G  A*  and  IF  =  Kr,  proving  the  following  corollary. 


r  G  A*. 

re  A*. 
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Corollary  5  The  irreducible  A-invariant  subspaces  of  KA  are  given  by  Kr, 
r  e  A *. 

Denote  the  complement  of  A  =  W  fl  A*  in  A*  by  Ac.  By  theorem  2  and 
theorem  4,  we  have  the  following  result. 

Corollary  6  KA  is  the  direct  sum  of  A-invariant  subspaces , 

KA  =  W®W\ 

where 

w  =  ®Kt = KAe'i 

rG  Ac 

and  e'  =  ±  £t€A<!  t. 

The  A-invariant  subspace  W  of  KA  in  theorem  4  can  also  be  described  by 

W  =  {/  €  KA  :  f(r)  =  0,  for  all  r  €  Ac}.  (1) 

In  general,  if  A  is  an  arbitrary  subset  of  A*  and  Ac  its  complement  in  A*,  the 
subspace  of  KA  define  in  ( 1)  is  A-invariant.  This  follows  from  the  diagonalizing 
formula  for  the  Fourier  transform  and  the  nonvanishing  of  characters. 

By  the  results  just  described,  single  characters  determine  irreducible  A- 
invariant  subspaces,  collections  of  characters  determine  A-invariant  subspaces 
and  partitions  of  the  collection  of  all  characters  determine  direct  sum  decompo¬ 
sitions  of  KA  into  A-invariant  subspaces. 

5.2  Nonsplit  case 

If  A  does  not  split  over  if,  then  KA  does  not  decompose  into  the  direct  sum  of 
one- dimensional  A-invariant  subspaces,  but  does  have  decompositions  into  the 
direct  sum  of  irreducible  A-invariant  subspaces.  For  K  —  Q  or  K  =  GF(pR)y 
p  a  prime  not  dividing  AT,  we  have  constructed  a  minimal  field  extension  of  K 
over  which  Cx  splits  or  equivalently,  having  a  primitive  N- th  root  of  unity.  If  N 
is  the  least  common  multiple  of  the  orders  of  the  elements  in  A,  then  the  same 
construction  produces  a  minimal  field  extension  E  of  K  over  which  A  splits. 
We  call  E  the  splitting  field  of  A  over  K . 

We  will  show  how  splitting  fields  of  A  over  K  can  be  used  to  construct 
decompositions  of  if  A  into  direct  sums  of  irreducible  A-invariant  subspaces. 
The  irreducible  A-invariant  spaces  will  not  generally  be  one-dimensional.  The 
first  new  concept  required  is  that  of  idempotents  which  replace  characters  as 
the  structural  basis  of  these  decompositions. 

Suppose  A  is  a  finite  abelian  group  of  order  N  and  if  is  a  field  over  which  A 
does  not  necessarily  split.  We  assume  throughout  that  if  =  Q  or  if  =  GF(pR) 
where  p  does  not  divide  N. 
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A  nonzero  element  e  €  KA  is  called  an  idempotent  if  e2  =  e.  Two  idempo- 
tents  ei  and  e 2  are  said  to  be  orthogonal  if  eie2  =  e2ei  =  0.  A  set  of  pairwise 
orthogonal  idempotents 

{ci,...,cj} 

is  called  a  complete  set  of  orthogonal  idempotents  if 

i=X> 

j=i 

The  sum  of  orthogonal  idempotents  e\  and  e2  is  an  idempotent  since 
(ei  4-  e*i)*  =  4*  eie2  4  62^1  4-  =  ci  4  62, 

and  if  e  is  an  idempotent,  the  set 

{e,l-e} 

is  a  complete  set  of  orthogonal  idempotents. 

Example  19  If  r  is  a  character  of  A  over  AT,  then  by  corollary  1, 


1 


is  an  idempotent  of  KA.  If  T\  and  r 2  are  characters  of  A  over  K,  then  the 
corresponding  idempotents  e\  =  and  e 2  =  ^^2  are  orthogonal. 

Example  20  Suppose  A  splits  over  K.  Denote  the  set  of  all  characters  of  A 
over  if-by  A*.  By  corollary  3,  the  set  of  idempotents 

{e  =  -ir:r€A*} 

is  a  complete  set  of  orthogonal  idempotents  of  KA. 

A  complete  set  of  orthogonal  idempotents  can  be  used  to  decompose  KA  into 
a  direct  sum  of  A-invariant  subspaces  in  much  the  same  way  as,  in  the  split  case, 
the  set  of  characters  produce  such  a  decomposition.  However,  the  A-invariant 
subspaces  generated  by  the  idempotents  will  not  generally  be  one-dimensional 
or  even  be  irreducible. 

Theorem  5  If 

(ej  :  1  <3<J} 

is  a  complete  set  of  orthogonal  idempotents,  then  KA  is  the  direct  sum  of  A- 
invariant  subspaces 

j 

KA  =  Y,®KAei- 

j= 1 
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Proof  Define 


Wj  =  KAeh  l<j<  J. 

is  the  A-invariant  subspace  of  KA  generated  by  tj.  Consider  the  sum 
£i=i  Wj.  By  orthogonality,  if 

j 

0  ”  ^  ^  fj  € 

i=i 

then 

0  =  ek  •  0  =  fk^k,  1  <  fc  <  J, 

and  the  sum  is  a  direct  sum  £j+1®Wj.  The  completeness  condition,  1  = 
£^=1  ej ,  implies  that  for  every  /  €  KA , 

j  j 

/  =  £/e;e£®Wi> 

i=i  i=i 

proving  the  theorem. 

The  ^-invariant  subspaces  If  Ae^,  1  <  j  <  J,  in  theorem  5  are  not  necessarily 
irreducible.  In  order  to  derive  conditions  on  and  eventually  construct  sets  of 
orthogonal  idempotents  producing  decompositions  into  irreducible  A-invariant 
subspaces,  we  will  use  the  splitting  field  E  of  A  over  K  and  the  Galois  group  of 
the  corresponding  extension  of  E  over  K. 

Suppose  that  E  is  the  splitting  field  for  A  over  K .  A  mapping  a  :  E  — ►  E 
is  callecTan  automorphism  of  E  over  K  is  a  is  an  automorphism  of  the  field 
E  fixing  the  elements  in  K.  The  collection  of  all  automorphisms  of  E  over  K 
forms  a  group  under  composition  called  the  Galois  group  of  E  over  K . 

Denote  by  T  the  Galois  group  of  E  over  K.  The  order  of  T  equals  the 
dimension  of  E  as  a  K-ve ctor  space. 

The  main  result  we  require  from  Galois  theory  is  the  following. 

Theorem  6  If  a  £  E  satisfies ,  for  all  a  £  F, 

a(a)  =  a, 

then  a  £  K. 

Suppose  that  c/>(N)  is  the  order  of  the  group  of  units  U(N)  of  Z/N . 

Example  21  The  cyclotomic  field  Q(£n),  =  e27ri^ ,  is  the  minimal  extension 

of  Q  over  which  Cn  splits.  Q(£n)  has  dimension  <j>{N)  as  a  vector  space  over 
Q.  The  Galois  group  of  the  extension  is  isomorphic  to  U(N)  having  elements 
uniquely  defined  by 

n  £  U(N). 
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Example  22  For  a  prime  p  not  dividing  N  and  R  the  smallest  positive  integer 
satisfying  pR  =  1  mod  N,  GF(pR)  is  the  splitting  field  of  Cn  over  GF(p). 
GF(pR)  has  dimension  R  as  a  vector  space  over  GF(p).  The  Galois  group  of 
the  extension  is  isomorphic  to  the  cyclic  group  Cr  and  is  generated  by  the 
automorphism  defined  by 

<j(a)  =  ap,  a  £  GF(pR). 


by 


T  acts  on  EA  and  on  =  ch(A,  E ).  For  /  £  EA  and  a  £  T,  define  fa  £  EA 

x  eA 


If  0*1,  cr2  €  then  fai<72  =  (Z^2)*1.  As  a  corollary  to  theorem  6,  we  have  the 
following. 

Corollary  7  For  f  £  EA,  fa  =  f,  for  all  a  £  F,  if  and  only  if  f  £  KA. 

For  r  £  A*,  ra  £  A*,  for  all  cr  £  F.  Define 


r t  =  {t°  :  a  £  r} 


and  call  Ft  the  T- orbit  at  r.  If  rf  £  Ft,  then  Ft  =  Ft'  and  if  IYi  ft  rr2  is  not 
empty  for  r2  £  A*,  then  Tri  =  Tr2.  The  action  of  F  on  A*  partitions  A* 
into  the  disjoint  union  of  T-orbits. 

A  subset  A  of  A*  is  called  F-invariant  if  r  £  A  implies  r°  £  A,  for  all  a  £  F. 
A  T-invariant  subset  A  of  A*  is  the  disjoint  union  of  the  T-orbits  contained  A. 

We  will  see  that  the  T-invariant  subsets  of  A*  determine  the  A-invariant  sub¬ 
spaces  of  KA  and  that  the  T-orbits  in  A*  determine  the  irreducible  A-invariant 
subspaces  of  KA.  Moreover,  we  will  show  the  relationship  between  decomposi¬ 
tions  of  T-in  variant  subsets  of  A*  into  disjoint  unions  of  T-invariant  subsets  and 
decompositions  of  A-invariant  subspaces  of  AT  A  into  direct  sums  of  A-invariant 
subspaces. 

For  r  £  A*,  define  eT  £  EA  by 


r'erv 


If  A  £  Ft,  then  eA  =  er  and  eT  depends  solely  on  the  T-orbit  at  r. 
Theorem  7  For  r  £  A *,  er  is  an  idempotent  in  KA. 


Proof  By  corollary  1,  e*  =  eT.  Since  r'  £  Fr  if  and  only  if  (r')*7  £  Tr  for  all 
cr  £  T,  we  have 

er  ~  eri  €  r, 

which  by  corollary  7  implies  er  £  KA  completing  the  proof. 
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Theorem  7  provides  a  method  for  constructing  idempotents  in  KA  from  T- 
orbits  of  characters  in  A*.  Moreover,  if  t*,  T2  €  A*  lie  on  distinct  T-orbits,  then 
by  corollary  1 

CTieT2  =  0) 

and  the  idempotents  en  and  eT2  are  orthogonal.  Each  T-orbit  determines  an 
idempotent  and  distinct  T-orbits  determine  distinct  orthogonal  idempotents. 

Denote  by  A* /Y  any  complete  set  of  representatives  of  T-orbits  in  A *.  A* /T 
is  a  set  formed  by  choosing  a  single  character  in  each  T-orbit.  If  A  is  a  T- 
invariant  subset  of  A*,  denote  by  A/T  the  subset  of  A* /T  formed  by  the  char¬ 
acters  contained  in  A. 

Theorem  8  The  set 

{eT  :  t  €  A* /T) 

is  a  complete  set  of  orthogonal  idempotents  in  KA. 

Proof  The  proof  follows  from  corollary  3  and 

reA'/V  reA* 


We  will  see  eventually  that  the  set  of  idempotents  in  theorem  8  is  special  in 
the  sense  that  the  A-invariant  subspace  generated  by  er,  KAeT)  is  an  irreducible 
A-invariant  subspace  in  KA. 

For  any  T-invariant  subset  A  of  A*,*define 


eA 


N  ^ 


re  A 


Since  A  is  T-invariant,  it  is  the  disjoint  union  of  T-orbits  and  we  can  write 
as  a  sum  of  orthogonal  idempotents 


^  ^  er, 

TGA/r 

proving  the  following  result. 

Theorem  9  For  any  Y -invariant  subset  A  of  A*, 


eA 


JL_ 

N 


Er 

re  A 


is  an  idempotent  in  K A  and 


KAe&  =  ^  @KAeT. 
reA/r 
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T-invariant  subsets  of  A  define  idempotents  in  KA  and  as  we  will  see  below 
every  idempotent  in  KA  comes  from  a  T-invariant  subset  by  the  formula  in 
theorem  9. 

Corollary  8  If  Ai  and  A2  are  nonintersecting  T-invariant  subsets  of  A* ,  then 
eAi  ond  eA2  are  orthogonal  idempotents  in  KA. 

Corollary  9  If  A  is  a  T-invariant  subset  and  Ac  is  its  complement  in  A* ,  then 
Ac  is  T-invariant  and 

{eA,eAc} 

is  a  complete  system  of  orthogonal  idempotents  in  KA. 

More  generally  if  A *  is  the  disjoint  union  of  T-invariant  subsets  A*  = 
Uf=1A a,  then  the  set 

{oAr  •  1  <r  <  R} 

is  a  complete  set  of  orthogonal  idempotents  of  KA.  By  theorem  5,  KA  is  the 
direct  sum  decomposition  of  the  ^-invariant  subspaces 

R 

KA  =  '£,®KAeAr. 

r—  1 

Corollary  10  For  T-invariant  subsets  Ai  and  A2  of  A*,  Ai  C  A2  if  and  only 
if  KAe^  C  KAe  a2. 

Proof  If  Ai  C  A2)  then  KAe Al  C  KAe&2  by  theorem  9.  Conversely,  if 
KAe&rJZ  KAe a2,  we  can  write 

OAi  fo a2j  f  ^  K A. 


For  any  r  G  Ai, 


r  =  reAi  =  fre  a2 


which  since  r  ^  0  implies  r  G  A2  completing  the  proof. 


Theorem  10  If  e  is  an  idempotent  in  KA,  there  exists  a  unique  T-invariant 
subset  A  of  A*  such  that  e  =  eA- 

Proof  Since  e  G  EA ,  we  can  write 

e = ^  E  ^r)T- 

t6>1* 

The  idempotent  condition  e2  =  e  implies 

e(r)2  =  e(r),  r  €  4*, 
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and  e(r)  is  either  0  or  1.  Denote  by  A  the  set  of  all  r  G  A*  such  that  e(r)  =  1. 
Since  ea  =  e,  for  all  a  G  I\  A  is  T-invariant  and  e  =  e^)  proving  the  theorem. 

Suppose  VT  is  an  A-invariant  subspace  of  if  A  and  WE  is  the  i?-subspace 
spanned  by  W  in  EA .  We  is  an  A-invariant  subspace  of  EA.  If 

{wr  :  1  <r<R} 

is  a  basis  of  the  if -vector  space  W ,  it  is  also  a  basis  of  the  jB-vector  space  We 
and  W  =  We  H  KA. 


Theorem  11  WE  is  F -invariant.  If  f  G  WE,  then  fa  G  WE)  for  all  a  G  T. 
Proof  If  /  G  WE)  then 

R 

f  =  5Z  «(r)  €  E, 

r=l 

and  for  all  cr  G  T, 

fi 

/*  =  a(r)au;r  G 

r—  1 

proving  the  theorem. 


Set  A  =  WeC\A*  and  denote  the  complement  of  A  in  A *  by  Ac.  Theorem  11 
implies  A  and  Ac  are  T-invariant.  By  theorem  4  especially  the  decomposition 
of  WE  given  in  (  1),  we  have 

W  =  {f€KA  :  f(r)  -  0,  for  all  r  G  Ac}.  (2) 

In  general  if  A  is  an  arbitrary  T-invariant  subset  of  A*  and  Ac  its  complement 
in  A*,  the  subspace  of  KA  defined  by  (  2)  is  A-invariant. 

Theorem  12  For  r  G  A*}  r  G  WE  if  and  only  if  er  G  W . 

Proof  Suppose  r  G  A*.  If  r  G  since  W#  is  T-invariant 

Tt  C  VTjs 

and  by  theorem  7  er  G  WE  O  KA  =  VT.  Conversely,  if  er  G  VT,  then  by 
corollary  1  and  the  A-invariance  of  WE 

1 

reT  =  —r  £We 

proving  the  theorem. 

Theorem  12  is  the  main  result  we  need  for  describing  A-invariant  subspaces 
of  if  A  in  terms  of  idempotents  of  if  A. 
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Corollary  11 


W  =  KAe A=  J2  ®KAer . 

rSA/r 


where  A  =  fl  A* . 

Proof  By  theorem  12,  eT  €  W ,  for  all  r  €  A,  implying  eA  €  W.  Since  W7  is 
A-invariant,  KAe a  C  W7.  Corollary  9  implies  1  =  eA  +  By  (  2),  for  any 

few, 

/  =  /•!  =  /eA  +  /eAc  =  /eA, 

completing  the  proof. 

By  corollary  11,  every  A-invariant  subspace  W  of  KA  is  generated  by  an 
idempotent  eA,  where  A  is  the  T-invariant  subset  We  H  A*.  By  corollary  10, 
we  have  the  following  result. 

Corollary  12  An  A-invariant  subspace  W  of  KA  is  irreducible  if  and  only  if 
W  —  KAer,  for  some  r  €  A*. 

The  preceding  discussion  reduces  the  study  of  A-invariant  subspaces  of  KA 
to  the  study  of  T-invariant  subsets  on  A*.  The  irreducible  A-invariant  subspaces 
of  KA  correspond  to  T-invariant  subsets  consisting  of  a  single  T-orbit.  A- 
invariant  subspaces  of  KA  which  are  the  direct  sums  of,  say  fc,  irreducible 
A-invariant  subspaces  correspond  to  T-invariant  subsets  consisting  of  exactly  k 
distineLT-orbits.  The  following  theorem  summarizes  these  remarks. 

Theorem  13  KA  is  the  direct  sum  of  irreducible  A-invariant  subspaces 

KA=  ®KAt. 

r<=A*/r 

The  number  of  irreducible  A-invariant  subspaces  in  K A  is  equal  to  the  number 
of  T -orbits  in  A* . 

Every  A-invariant  subspace  W  of  K A  is  the  direct  sum  of  irreducible  A- 
invariant  subspaces 

W=  Y,  ®KAt, 
r€A/r 

where  the  number  of  irreducible  A-invariant  subspaces  in  W  is  equal  to  the 
number  ofT-orbits  in  We • 
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Suppose  G  is  an  arbitrary  finite  group  of  order  N  and  if  is  a  field.  Usually 
in  examples  and  applications,  if  =  Q,  R,  C  or  a  finite  field  GF(pR ),  p  a  prime 
not  dividing  N,  but  results  hold  generally  for  fields  whose  characteristic  does 
not  divide  N, 

The  group  algebra  KG  of  G  over  K  is  the  if- vector  space  of  all  formal  sums 

/  =  £/(*)*,  meK, 

t£G 

under 

/  +  0  =  £(/(*) +  $(*))*,  /,  s  e  KG, 

t^-G 

a/  =  E(a-fW)f) 

with  if -algebra  multiplication 

fg  =  E  (e  f(u)g(u~ 

t£G  Vti^G 

The  additions  and  multiplications  inside  the  brackets  are  taken  in  if. 

Generally  since  G  is  not  necessarily  abelian,  fg  is  not  necessarily  equal  to 
gf,  f)  9  €  KG.  Replacing  u  by  tu~l  in  the  inner  summation, 

/0  =  E  ( E  0(u)/(*u_1) 

— t£G  \u6G 

Denote  by  L{G\  if)  the  if -vector  space  of  all  if-valued  functions  on  G.  Every 
/  G  L(G ;  if)  defines  a  formal  sum  in  KG 

/  =  £/(*)* 
t6G 

and  we  can  identify  the  if-vector  space  KG  with  the  if-vector  space  L(G ;  if). 
The  multiplication  fg  in  KG  corresponds  to  the  possibly  noncommutative  con¬ 
volution  in  L(G\  if) 

/  *  9(t)  =  E  /(w)0(«_1O  =  E  5(u)/(^_1).  /.  0  €  L(G;  X). 

Under  this  identification  the  delta  function  <5W,  u  £  G,  corresponds  to  a  formal 
sum  having  a  single  nonzero  coefficient  which  we  denote  by 

u  =  Y6u(t)t 

teG 


/,  0  e  KG. 


f  €  KG,  a  GK, 


t)  It,  f,ge  KG. 
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In  this  way  we  can  view  G  as  a  subset  of  KG  and  a  basis  of  the  K- vector  space 
KG. 

For  u  (E  G  and  /  €  KG) 

«/ = Yj  f^ut  =  /(u_10f- 

t£G  t£G 

We  call  uf  the  left  translation  of  /  by  u.  Right  translation  can  also  be  defined 
and  generally  differs  from  left  translation,  unless  G  is  abelian.  For  u  €  G,  the 
operator  L(u)  of  KG  defined  by 

L(u)f  =  uf,  fe  KG , 

is  a  linear  isomorphism  of  the  K- vector  space  KG  and  we  have 
L(uv)  =  L(u)L(v ),  u,  dGG, 

L(w_1)  =  L(u)-1,  «  €  G, 

where  L(u)L(v)  and  L(u)~1  denote  composition  and  inversion  of  operators  on 
KG. 

Denote  by  GL(KG ,  K)  the  group  of  linear  isomorphisms  of  the  K-ve ctor 
space  KG.  The  mapping 


L  :  G  — ►  GL(KG,K) 

is  an  isomorphism  of  the  group  G  into  the  group  GL(KG,K )  called  the  left- 
regular  representation. 

A  subspace  V  of  the  if -vector  space  KG  is  called  G- invariant  if  for  all  u  €  G, 
uV  =  {uf:f€V}c  V ; 

A  G-invariant  subspace  V  of  KG  is  called  irreducible  if  the  only  G-invariant 
subspaces  of  V  are  (0)  and  V.  One  of  the  main  goals  of  nonabelian  group  har¬ 
monic  analysis  is  to  characterize  the  G-invariant,  irreducible  G-invariant  and 
direct  sum  decompositions  of  KG  into  irreducible  G-invariant  subspaces.  For 
an  abelian  group,  its  character  theory  provided  all  the  necessary  tools  for  an¬ 
swering  these  questions.  Generally  these  problems  require  a  vast  mathematical 
machinery  for  their  solution.  However  for  the  nonabelian  groups  considered  in 
this  text,  explicit  solutions  for  many  of  these  questions  will  be  derived  in  terms 
of  abelian  group  character  theory. 

For  g  €  KG ,  the  operator  L(g)  of  KG  defined  by 

L(g)f  =  gf,  f  e  KG, 

is  a  linear  homomorphism  of  the  if-vector  space  KG.  Since 
L(g)  =  '£/g(t)L(t)>  g  €  KG, 

t£G 
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G- invariant  subspace  V  of  KG  are  if  G-invariant  subspaces.  If  V  is  a  G-invariant 
subspace  of  KG  and  /  £  V,  then  for  all  g  £  KG, 

gV  c  v. 

Idempotent  theory  provides  a  convenient  language  in  which  to  express  many 
of  the  results  in  the  following  chapters.  However  explicit  results  can  usually  be 
written  in  terms  of  abelian  group  character  concepts. 

A  nonzero  element  e  £  KG  is  called  an  idempotent  if  e2  =  e.  Two  idempo- 
tents  e\  and  in  KG  are  called  orthogonal  if  e^  =  e2ei  =  0.  A  set  of  pairwise 
orthogonal  idempotents 

{ej  :1<3<J} 

is  said  to  be  complete  if 

i  «5> 

i 

Example  23  If  e  is  an  idempotent  in  KG,  then  {e,  1  —  e}  is  a  complete  set  of 
orthogonal  idempotents. 

Example  24  If  e  is  an  idempotent,  then  the  G-invariant  subspace  generated 
by  e,  KGe,  has  e  as  a  right  unit 

KGe  =  {aG  KG  :  ae  =  a}. 

Simply  write  a  £  KGe  as  a  =  a'e,  a'  €  KG  and  use 

ae  =  a'e2  =  a'e  =  a. 

Conversely  if  a  £  KG  can  be  written  as  a  =  ae,  then  ae  £  KGe  implies 
a  £  KGe. 

Example  25  If  e  is  an  idempotent  then  KG  =  KGe  ®  KG(  1  —  e).  Since 
1  =  e  +  (1  —  e),  for  any  a  £  KG 

a  =  a  •  1  =  ae  4-  a(l  —  e)  £  KGe  4-  K G(1  -  e) 

and  KG  =  KGe  +  KG(  1  —  e).  The  sum  is  a  direct  sum  since  if  a  £  KGe  fl 
KG(  1  —  e)  and  we  write  a  =  a'e  =  a"(l  -  e),  then 

a  =  a'e  =  a"(l  —  e)e  =  a"(e  —  e2)  —  0. 


The  result  described  in  example  25  holds  generally  for  any  complete  set  of 
orthogonal  idempotents. 
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Theorem  14  If 


{ej  :  1  <j<J} 

-  is  a  complete  set  of  orthogonal  idempotents,  then 

j 

KG  =  '±2®KGej. 

3- 1 

Proof  Since  1  =  ej>  f°r  any  a  €  KG, 

a  =  Y.aei  £j2KGei 

3=1  3=1 

and  KG  =  KGe,.  If 

J 

0  =  iej9 

i=i 

then  for  any  1  <  k  <  J,  by  orthogonality 

j 

0  =  0  •  ek  =  ]T  OL^ejGk)  =  akek 
j= 1 

proving  the  sum  is  a  direct  sum.  The  same  argument  shows  the  following  corol¬ 
lary.  _ 

Corollary  13  If  an  idempotent  e  can  be  written  as  the  sum  of  two  orthogonal 
idempotents,  e  =  e\  +  e2,  then 

KGe  =  KGex  ©  KGe2. 


Example  26  If  KG  =  W\  0  W2j  where  W\  and  W2  are  G-invariant  subspaces 
and 

1  =  Ci  4*  ^2,  e\  G  W\ ,  e2  G  W2, 
then  {ei,e2}  is  a  complete  set  of  orthogonal  idempotents.  Since 

ei  =  ei  •  1  =  e\  +  e\e2 

with  ei,  e\  €  W\  and  eie2  G  W2,  by  uniqueness  of  representation  e\  =  e\  and 
eie2  =  0.  Now 

ei  =  1  •  ex  =  e?  +  e2ei  =  ci  +  e2ei 
implies  ^ei  =0  verifying  the  example. 
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Generally  we  have  the  following  result. 

Theorem  15  If  KG  is  the  direct  sum  of  G -invariant  subspaces 

J 

KG  =  Y/®Wj 

j  =  l 

and 

j 

1  =  53  eh  ei  €  wi>  1  <j<J, 

i=i 

then 

{ej  :  1  <  j  <  J} 

is  a  complete  set  of  orthogonal  idempotents. 

Proof  Completeness  is  by  definition.  For  any  1  <  k  <  J, 

J 

eh  —  6^  •  1  =  ^  ^ 
j= i 

where  e^ej  €  Wj,  1  <  j  <  J}  which  by  uniqueness  of  representation  implies 

e*  =  e\,  ekej  =  0 ,  j^k 

competing  the  proof. 


Example  27  If  KGe  =  W\  ©  W2,  the  direct  sum  of  G-invariant  subspaces  W\ 
and  W2  with  e  an  idempotent  and  e  =  ex  +  e2,  ex  €  WX)  e2  €  W2,  then  ex  and 
e2  are  orthogonal  idempotents.  By  example  24 

"  ei  —  eie  =  e\  -h  eie2 

which  by  uniqueness  of  representation  implies 

ei  =  ef,  eie2  =  0. 

The  same  argument  shows  that 

e\  =  e2,  e2e\  =  0. 


Generally  we  have  the  following  which  we  state  without  proof. 
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Theorem  16  If  KGe  is  the  direct  sum  of  G -invariant  subspaces 

j 

KGe  =  '%2®Wj) 

j- 1 

with  e  an  idempotent  and 


J 

e  ~  5Z  ei  ^  ^  ’  *  -  J  -  J) 

j= i 

then 

{ej  :l<j<J} 

is  a  set  of  pairwise  orthogonal  idempotents. 


A  homomorphism  P  of  the  If -vector  space  KG  is  called  a  projection  if  P 2  = 
P.  If 

imP^{P(a):aelfG} 
kerP  =  {a  €  ifG  :  P(a)  =  0}, 
then  imP  and  kerP  are  subspaces  of  KG  satisfying 


im  P  =  {a  e  KG  :  Pa  =  a} 

KG  =  im  P  ©  ker  P. 

Every  subspace  W  of  if  G  determines  a  projection  P  such  that  W  =  im  P.  Such 
a  projection  exists  since  if  ifG  =  W  ®  W',  for  some  subspace  W*  of  KG ,  then 
we  can  define  P  by 


P(u>  -f  u/)  =  w,  w  €  W,  it;'  €  W7. 

We  will  now  show  that  every  G-invariant  subspace  W  of  K G  is  generated 
by  an  idempotent.  Consider  any  projection  P  of  ifG  satisfying  W  =  imP  and 
define  the  mapping  P0  :  IfG  — *  ifG  by 

Po(a)  =  jf  Y,  u_lp(u°)>  «  € 

Po  is  a  homomorphism  of  the  if-vector  space  KG.  Since  P(ua)  €  W,  «EG 
and  a  €  KG,  and  W  is  G-invariant 

P0{a)eW,  a  £  KG. 
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Theorem  17  For  all  a  €  KG 


and  if  a  £W,  then 


P0(a)=aP0(l), 


Po(a)  =  a. 


Proof  For  t  E  G, 

t-'Poit)  =  t-'v-'PM  =  Tf  £  u"lp(u)  =  p0(l), 

UGG  ti€G 

by  a  change  of  variables.  The  linearity  of  Po  proves  the  first  part.  If  a  E  W, 
then  ua  E  W,  u  E  G  and  P(ua)  =  ua  implying 

Po(ua)  —  ^  u~lua  =  a, 

N  ueG 

completing  the  proof. 

Since  im  Po  C  W  and  Pq  acts  by  the  identity  mapping  on  W,  by  theorem  17 
we  have  the  following. 

Corollary  14  Po  is  a  projection  of  KG  satisfying  W  ~  imPo- 
SetT^=  Po(l). 

Corollary  15  e  is  a  generating  idempotent  for  W, 

W  =  KGe ,  e2  =  e. 


Proof  By  theorem  17  and  corollary  14,  for  all  w  E  W> 

w  =  P0{w)  =  we . 

Since  e  E  W,  e  is  an  idempotent  in  W  and  a  right  unit  for  W.  Example  24 
implies  W  =  KGe ,  completing  the  proof. 

Since  every  G-invariant  subspace  W  has  a  generating  idempotent,  by  corol¬ 
lary  15,  there  exists  a  G-invariant  subspace  Wf  such  that  KG  ~  W  ©  W'.  More 
generally  we  have  the  following  result. 
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Corollary  16  If  W\  and  W 2  are  G-invariant  subspaces  of  KG  such  that  W\  C 
W2  then  there  exists  a  G -invariant  subspace  W[  ofW2  such  that  W2  =  W\  ®  W[ . 


Proof  If  Wf  is  a  G-invariant  subspace  of  KG  such  that  KG  =  W\  ©  W' ,  then 
W2  =  W\  ©  (W'  fl  W2)  where  W{  =  Wf  ft  W2  is  a  G-invariant  subspace  of  W2. 

An  idempotent  e  E  KG  is  called  primitive  if  e  can  not  be  written  as  the  sum 
of  orthogonal  idempotents  in  KG. 

Theorem  18  e  is  a  primitive  idempotent  in  KG  if  and  only  if  KGe  is  irre¬ 
ducible. 

Proof  If  e  is  not  primitive  and  e  —  e\  -f  e2,  where  e\  and  e2  are  orthogonal 
idempotents,  then  KG  =  ifGei0ifGe2  and  KGe  is  not  irreducible.  Conversely 
if  W  is  a  G-invariant  subspace  of  KGe,  then  by  corollary  16,  KG  —  W®  W',  for 
some  G-invariant  subspace  Wf  of  KGe.  Example  27  implies  e  is  not  a  primitive 
idempotent,  completing  the  proof 

The  problem  of  constructing  G-invariant,  irreducible  G-invariant  and  di¬ 
rect  sum  decompositions  of  KG  into  irreducible  G-invariant  subspaces  can  be 
replaced  by  that  of  constructing  idempotents,  primitive  idempotents  and  com¬ 
plete  sets  of  primitive  orthogonal  idempotents. 

Decompositions  of  KG  into  direct  sums  of  G-invariant  subspaces 

J 

KG  =  '£j@Wj 

—  j= 1 

lead  to  block  diagonal  matrix  representations  of  the  left  translation  operators 
L(a),  a  e  KG ,  and  to  fast  algorithms  for  computing  products  in  KG.  Details 
for  special  groups  G  will  be  given  in  the  following  chapters.  However,  in  this 
chapter  we  will  describe  the  general  outline  of  the  method. 

For  t  £  G,  the  matrix  S(t)  of  L(t)  relative  to  the  basis  G,  of  KG ,  ordered 
in  some  way,  is  a  permutation  matrix  reflecting  the  group  structure  of  G.  The 
matrix  5(a)  of  L(a),  a  E  KG ,  relative  to  G  is  a  linear  combination  over  K  of  the 
permutation  matrices  S(t ),  t  E  G.  If  G  is  a  cyclic  group,  then  the  corresponding- 
mat  rices  5(£),  t  E  G,  are  cyclic  shift  matrices  and  the  corresponding  matrices 
S(a),a  €  KG,  are  circulant  matrices.  In  this  case  the  direct  sum  decomposition 
of  KG  given  by  the  characters  diagonalizes  the  circulant  matrices  with  the 
Fourier  transform  describing  the  change  of  basis. 

Generally  since  Wj  is  G-invariant,  L(a),  a  £  KG,  maps  Wj  into  itself. 
Denoting  by  Tj(a)  the  matrix  of  the  restriction  of  L(a)  to  some  basis  of  Wj, 
the  matrix  direct  sum 

T(a)  =  ^rj(a), 

3= 1 
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is  the  matrix  representation  of  L(a)  relative  to  the  basis  of  KG  formed  by 
piecing  together  the  bases  from  the  Wj,  1  <  j  <  J.  We  say  the  new  basis  is 
compatible  with  the  direct  sum  decomposition  of  KG. 

To  compute 

a/3  =  L(a)/3 ,  a,  (3  e  KG , 

we  write  /3  in  terms  of  the  compatible  basis,  compute  T(o)  and  form  the  matrix 
product  of  T(a)  with  the  coordinates  of  /?  in  the  compatible  basis  and  then 
translate  the  result  back  to  the  basis  G. 
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6  Group  algebras  of  A<)  B 

A  subgroup  A  of  a  group  G  is  called  normal  if  for  all  t  G  G, 
tAt~l  =  {tat~l  :  a  G  A}  C  A. 

G  is  said  to  be  the  semidirect  product  of  a  normal  group  A  and  a  subgroup  B  if 
every  t  e  G  can  be  written  uniquely  as 


t  =  xy ,  x  G  A,  y  G  B. 


In  this  case  we  write  G  =  A<$B  and  we  call  x  the  A- component  and  y  the 
B-component  of  t.  If  t  =  xy  and  t(  =  xV  with  x,  x'  e  A,  y>  y'  e  B  then 

ttf  =  x{yx,y^l)yyt 

where  x(yxfy~l)  is  the  A-component  and  yyf  is  the  B  component  of  ttf , 

In  this  chapter  we  will  study  groups  of  the  form  G  =  A<£  B  where  A  and  B 
are  abelian  groups. 

For  a  G  KA  and  /3  G  we  can  view  a,  (3  €  IfG.  If  a/3  =  0,  then 


2  H  ot(x)/3{y)xy  =  0 

xeAyeB 


which  implies 
and  a  =  0  or  /3  =  0. 


a(x)/3(y)  =  0,  x  e  A>  y  €  B, 


6.1  Split  case 

Suppose  G  =  A  B  where  A  and  B  are  abelian  groups  splitting  over  K.  Denote 
the  character  groups  of  A  and  B  over  K  by  A*  and  B *.  If  r  €  A*  and  A  G  B*, 
then  as  elements  in  KG 

t  =  E  t{x)x,  t(x)  €  K, 
xeA 

and 

a  =  E  A(y)y>  e  K- 

yeB 

B  acts  on  A*.  For  r  G  A*  and  y  G  B  define  ry  G  KG  by 

Ty  =  yry~1. 

Since  A  is  a  normal  subgroup  of  G 

ry  =  E  T(x)yxy~'1  =  E  T{v~lxy)x 

x£A  xEA 
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is  also  in  A *.  For  y,  y'  £  B 


(ry)y/  =  ry,y 

and  as  a  special  case 

( Ty)y _1  =  r. 

For  r  £  A* 

Bt  =  {ry  :  y  £  B} 

is  a  subset  of  A *  called  the  J3-orbit  over  r. 

If  r  and  r'  are  in  -4*  and  r'  €  .Br  with  r'  =  ry,  */  G  £,  then 

r  =  (ry)y  1  =  r/y  €  Brf 

and  Br  =  .Br'. 

For  r  and  r'  in  A*,  if  £?r  fl  Br '  ^  with  ry  ~  r/y  ,  y,  y'  G  then 
r  =  (ry)y  1  =  (^(r')y/ j  =  rtV  y 

and  r  €  Br*  implying  Br  =  Br'.  Consequently  A*  is  partitioned  into  the 
disjoint  union  of  distinct  B-orbits.  A  set 


C  A* 


is  called  a  complete  set  of  representatives  for  the  collection  A* /B  of  distinct 
j3-orbits  in  A*  if  A*  is  the  disjoint  union 


A* 


K 


Y,Brk- 

k=  1 


Since 


‘■{E 


T, 


1  = 


M 


E* 

AeB* 


where  L  and  M  are  the  orders  of  .4  and  B 


N 


E  E>*. 


where  AT  —  LM  is  the  order  of  G .  However  generally  r  and  A  do  not  commute. 
We  do  have  the  following. 


Theorem  19  For  t  =  xy  £  G,  x  £  A,  y  £  B, 

*=  jf  ]£  X)  TS'(^_1)%_1)Ty^ 

V  r6A*  A€B* 
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Proof  Since  xry  =  ry(x~l)ry  and  yX  =  X{y~1)X) 

trX  =  xyrX  —  xryyX  =  ry  (x^1)X(y~l)ry  X 
completing  the  proof. 

Since  G  is  a  basis  of  the  K-ve ctor  space  KG ,  by  theorem  19,  the  collection 
of  products  in  KG 

{tX'.t£A\  X  G  B*} 

is  also  a  basis  of  KG.  Generally  KGrX  is  not  one-dimensional,  so  that  the 
spaces 

{KGrX  :  r  G  A*9  X  G  B*} 

intersect.  If  G  is  abelian,  then  the  products  rA,  t  G  A*,  A  G  £*,  are  the 
characters  of  G  over  If  and  the  results  of  chapter  4  apply. 

Generally  to  form  direct  sum  decompositions  of  KG  into  left  G-invariant 
subspaces  we  must  modify  the  above  approach.  For  r  G  A*,  define 

B{t)  =  {y  €  S  :  rv  =  r}. 

B(r)  is  a  subgroup  of  B  called  the  centralizer  of  r  in  B. 

Suppose  Ti,  t2  €  A*  with  r2  =  rf1,  yx  €  B.  If  y  €  B(rx ),  then  since  B  is 
abelian 

Ty  =  (Tyi)y  =  (ry)yi  =ryi  =r2 

and  y  €  £(t2)  proving  the  following  result. 

Theorem  20  J/r2  G  t1?  t2  G  A*,  then 

B{n)  =  b[t2). 

We  can  assign  a  centralizer  in  B  to  every  B-orbit  in  A*.  The  assumption 
that  B  is  abelian  is  essential. 

Suppose  r  G  i4*  and  consider  B(r).  For  y  G  B)  the  set 

yB(r)  =  {yz  :  z  G  B(r)}  C  B 

is  called  the  left  coset  of  B(r)  in  £  determined  by  y.  If  y'  G  £(r),  then 
y'£(r)  =  y£(r).  The  usually  arguments  show  that  two  left  cosets  are  either 
equal  or  else  have  empty  intersection.  The  collection  B/B(r)  of  left  cosets  of 
B(r)  in  B  forms  a  partition  of  B.  A  set 

[ys  :  1  <  s  <  S}  C  B 

is  called  a  complete  system  of  representatives  of  B/B{r)  if  B  is  the  disjoint 
union 

s 

B  =  X>S(r). 

5=1 
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Since  B  is  an  abelian  group,  a  group  structure  can  be  placed  on  B/B(t)  by 

yB{T)y'B(T)  =  yy'B(r ) 

{yB{T))-ly~l B{t),  y,  y'  €  B. 

For  y  G  y5B(r),  with  y  =  ysz ,  1  <  s  <  S,  z  G  B(r), 

Ty  =  (r*)v'  =ry%  1  <  5  <  5. 

Moreover  if  rVa  =rVt,  1  <  s,  £  <  S,  then 

rvrlVs  _  T 


and  ysB(r)  =  ytB(r).  Consequently  the  B-orbit  Br  has  order  S  and 

Br  =  {rVs  :  1  <  s  <  S'}. 

Generally  the  ys,  1  <  s  <  S  depend  on  r  and  when  we  need  to  express  the 
dependence  we  write  yj,  1  <  s  <  ST.  If  Mr  denotes  the  order  of  B(r), then 
M  =  MrSr,  where  M  is  the  order  of  B. 

We  will  now  show  that  the  collection  of  products 

~rX,  re^AeB(T)*, 

forms  a  complete  set  of  primitive  orthogonal  idempotents  for  KG.  We  break  up 
the  proof  into  the  following  three  theorems. 

Theorem  21  For  Ti,  t2  G  A*  and  Ax  G  B(rx)*,  A2  G  B(t2)*, 

(TiAi)2  =LMnr1Ai 

and 

(*ri  Ai  )(t*2 A2)  =  0,  unless  rx  =  r2  and  Ax  =  A2. 


Proof  Since  rx  =  t1}  y  €  B(r\), 

T1A1  =  ^  Ai(3/)'riy  =  ^  ^i(s/)j/Ti  =  A1T1. 

y€B(n)  y€B(n) 

Consequently 

(nAi)2  =  Ait2Ax  =  LtiA2  =  LMn  7*iAi 

and 

(riAi)(r2A2)  =  Ai(tit2)A2  =  0,  unless  n  =  r2, 
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in  which  case 


(tiAi)(t1A2)  =  T1A1A2  =  0,  unless  Ai  =  A2, 
completing  the  proof. 

By  theorem  21,  the  collection  of  products  is  a  set  of  orthogonal  idempotents  in 

KG. 

Theorem  22 


1  =  Y  T 

'  ki->.ktrYLM' 


Proof  The  theorem  follows  from 


‘■iE 


and 


tga * 


Theorem  22  implies  completeness. 

For  r  £  A*  and  A  £  B(r)*,  by  theorem  21 

-  e-mTX 

is  an  idempotent.  Since  j^T  is  an  idempotent  and  r\  —  At 

1  1 

e  =  —re  =  —  er. 

1/  L 

Theorem  23  If  r  £  A*  and  A  £  then  e  =  is  a  primitive  idem- 

potent. 


Proof  Assume  e  is  not  primitive  and  e  =  ei+e2  where  e\  and  e2  are  orthogonal 
idempotents.  Since 

1  1 


and 


we  have 


ex  =  e\e  =  — e^r  =  —e\T 
1/  1/ 


d  =  eei  =  yree  1  =  —rei 
1/  I 


2  1  2  1 

el  =  el  =  ^2relr  =  prelT- 
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ei(r',2/)  €  K. 


Since  A*  is  a  basis  of  K A ,  we  can  write 


J2  (  ei(r',3/)T/N)  y, 

y£B  Vr'GA-  / 


Consequently 

However 

and 

implying 


ei  =  j?TeiT  =  73  H  H  e1(r',y)Tr,j/r. 


L3 


y^Br'eA* 


rr'yr  =  0,  unless  r'  =  r  and  $/  E  H(r), 

T2yr  =  L2Ty,  y  €  B(r), 

ei  =  yr  ei(r,y)y  =  yrmi,  mi  €  KB(t). 


V€B(t) 

The  same  argument  shows  that 


e2  =  yrm2,  m2  €  KB(t), 

Ll 


implying 

e=is;TA=r<m'+m2)' 

Since  flT  KA  and  A,  mi  -f  m2  E  O(r), 


——A  =  mi  -f  m2. 
Mr 


We  will  show  that  m\  and  m2  are  orthogonal  idempotents  in  KB(t),  contra- 
dieting  the  fact  that  jj-  A  is  a  primitive  idempotent  in  KB(r). 

Since  rmi  =  mix, 


1 

TTm  1 


implying  =  m\*  The  same  argument  shows  m\  =  m2  and  mim2  =  m2mi  = 
0,  completing  the  proof. 


We  have  proved  that  the  collection  of  products 

jl-rA:r^-,AcB(rr 

is  a  complete  set  of  primitive  orthogonal  idempotents  for  KG  and  we  can  apply 
the  results  of  the  previous  chapter. 
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Theorem  24 


Y  (&KGt\ 

\£B(Ty 

with  the  direct  sum  factors 

KGrX ,  tgA\  A  €  B{t)\ 
irreducible  G-invariant  subspaces  of  KG . 

In  the  abelian  case,  the  irreducible  G-invariant  subspaces  of  KG  are  one¬ 
dimensional  and  are  uniquely  determined  by  the  characters.  The  Fourier  trans¬ 
form  is  the  change  of  basis  transform  between  the  G-basis  and  the  character 
basis.  The  importance  of  the  Fourier  transform  in  digital  signal  processing  is 
greatly  enhanced  by  fast  algorithms  for  the  Fourier  transform. 

In  contrast,  if  G  is  not  abelian,  an  irreducible  G-invariant  subspace  of  KG 
is  not  necessarily  one-dimensional.  For  groups  of  the  form  G  =  A<$B,  with 
A  and  B  abelian,  we  will  determine  bases  for  the  subspaces  KGrX ,  r  £  A*) 
X  £  B(r )*,  and  derive  fast  algorithms  for  relating  components  in  the  G-basis  to 
components  in  the  new  basis. 

Suppose  r  £  A*  and  A  £  B(r )*.  Set  y8  =  yj,  1  <  s  <  S  =  ST  in  the 
following  discussion.  For  x  £  A  and  y  £  £?,  with  y  =  y3z ,  1  <  s  <  S  and 
z  e  b(t), 


KG  =  Y, 


— ~  xyrX  —  xyszrX  =  xrVayszX 

=  rVa  (a:”1 )  A  (z"1  )ysrX. 

Theorem  25  For  r  £  A*  and  X  £  B(r)*,  the  set 

{ysr X  :  1  <  s  <  S} 

is  a  basis  of  the  K -vector  space  KGrX. 

Proof  Since  the  set 

{xyrX  :  x  £  A,  y  £  B} 


spans  KGrX ,  the  set 


{ysrX  :  1  <  s  <  S} 

spans  KGrX.  We  must  show  that  this  set  is  linearly  independent.  Suppose  that 
5 

0-E a(s)y$rX,  a(s)  £  K,  1  <  s  <  S. 

5=1 
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Multiplying  on  the  left  by  rVt ,  1  <  t  <  S', 

s 

0  =  rf  0  =  ^2  ot{s)rytrVoys  A 

5=1 

=  La(t)ytT  A 

implying  a(<)  =  0.  Since  £  is  arbitrary,  1  <  £  <  S,  a(t)  =  0,  for  all  1  <  £  <  5, 
completing  the  proof. 


We  will  now  use  theorem  26  to  derive  an  algorithm  for  computing  the  coef¬ 
ficients 

aT\(s),  t  €  A*,  A  6  B(t)*,  1  <  s  <  ST, 

based  on  abelian  group  Fourier  transforms. 

For  y  e  B,  define  ay  €  KA  by 


Qy(x)  =  a(xy ),  x  €  A, 
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and  compute  the  Fourier  transform  of  ay  over  -4, 

«»(T)  =  J  Y  “yfcM®-1)*  T  6  A*- 

M  Fourier  transforms  over  A  are  required  at  this  stage.  Place  the  results  into 
a  two-dimensional  array  over  A*  x  B 

a(r,  y)  =  a„(r),  rZA*,y£B. 

For  each  r  G  j4*  we  will  compute  the  coefficients 

ctrx(s),  1<s<St,\€ 

For  1  <  s  <  ST,  define  7J  €  KB(t)  by 

is  (z)  =  7(T,/S .  y*z)>  z  e  B(r)»  Vs  =  yl, 

and  compute  the  Fourier  transform  of  7J  over  B(r) 

7j(A )  =  -±r  Y  7:(*)A(*-1),  A  e  B(t)*. 

r  z€B(t) 

Sr  Fourier  transforms  over  B(r)  are  required  in  this  stage.  Implementing  this 
stage  as  r  runs  over  A *  requires  $ T  Fourier  transforms  over  abelian 

groups  over  varying  sizes. 

By  theorem  26,  for  each  r  £  A*, 

Q:ta(s)  —  7j(^)»  ^  >  1  ^  ^  <  5r. 


First  Stage  a  £  KG 

OLy  £  KA ,  2/  €  B 

FT  over  A,  y  G  B. 
at  £  Xj4*,  2/  G  B 

d  €  L(i4*  x  B) 
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Second  Stage  r  €  A 


d  €  L{A *  x  B) 

7J  £  KB(t),  1  <  s  <  ST 

FT  over  B{t),  1  <  s  <  ST. 

7j  €  KB(t)*,  1  <s<ST 


otr\(s)  =  71(A),  A  B(t)*,  1  <s<ST. 


The  algorithm  computes  the  coefficients  ; 


ar*(s),  r  e  A*,  A  €  1  <  s  <  Sr 

by  firsfTomputing  M  Fourier  transforms  over  A  and  then  partitioning  the  re¬ 
maining  computation  into  L  parallel  stages  parameterized  by  r  €  i*.  For  each 
r  e  A*,  we  implement  a  data  rearrangement  step  which  forms  ST  elements  in 
KB(t)  and  then  computes  ST  Fourier  transforms  over  B(r).  For  A  e  B(t)*> 
the  Sr-coefficients 

1  s  ^ 


are  given  by  evaluation  of  the  ST  Fourier  transforms  at  A. 

Assuming  that  the  Fourier  transform  over  an  abelian  group  of  size  P  requires 
PlogP  complex  add-multiplies,  the  first  stage  requires  N  log  L  complex  add- 
multiplies  to  implement  M  Fourier  transforms  over  A.  For  each  r  €  A*,  the 
second  stage  requires  MlogMr  complex  add-multiplies,  since  M  =  SrMT}  to 
implement  ST  Fourier  transforms  over  B(t ).  To  implement  the  second  stage 
over  all  r  £  A*  requires  M  J2T€Am  1°6  complex  add-multiplies.  The  complete 
algorithm  requires 

NlogL  +  M  ]T  log MT 

t£A* 


complex  add-multiplies. 
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Up  to  scale  multiple,  the  abelian  Fourier  transform  can  be  realized  by  a 
unitary  matrix  which  relates  the  coefficients  over  G  with  the  coefficients  over 
the  character  basis.  Since  the  relationship  between  the  coefficients  over  G  with 
the  coefficients  over  the  basis 

{s/JrA  :  r  €  A*,  A  €  B(t)*  ,  1  <  s  <  ST} 

is  based  on  abelian  group  Fourier  transforms  we  can  represent  the  mapping 
between  these  coefficients  by  a  unitary  matrix,  up  to  scale  multiple,  as  well. 
The  algorithm  can  be  displayed  by  the  following  four  pictures. 
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Data  readdressing 


aar(z)  =  a(ry’,ysz),  z  G  B(t),  1  <s<ST- 


VsB{t) 


2/2  B(t) 


7? 


yi  b(t) 


S  to 


Fourier  transform  over  B(r). 


7t 


B(t) 


'll 


B*{t) 


B{r) 


B*{r) 


7r 


B(t) 


B-(t) 


Data  readdressing 


H - b*(t) 

X 

4 -  B*(t) 

A 


- 1 -  B*(t) 

X 

aTX(s)  =  7?(A),  l<S<Sr,reA*,Ae  B{r)\ 
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From  the  first  data  rearrangement  stage  we  can  see  how  non-commutativity 
affects  the  classical  abelian  group  two-dimensional  Fourier  transform  A  x  B. 
If  B(t)  =  B,  then  ST  =  1  and  the  next  stage  Fourier  transform  is  over  all 
of  B.  For  B(t)  =  (0),  S  =  M)  the  order  of  B)  and  the  data  passes  through 
unchanged.  No  Fourier  transform  is  required.  If  B(t)  is  a  proper  subgroup  of  B, 
the  Fourier  transform  of  the  decimated  data  over  B(r)  produces  a  periodization 
of  the  Fourier  transform  over  B .  Generally  the  coefficients  aT\(s)  represent 
periodizations  and  phase  modulated  periodizations  of  the  Fourier  transform 
over  B  of  the  Fourier  transform  of  the  initial  data  with  respect  to  A. 
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7  Examples 

7.1  The  Dihedral  Group  D ^  =  CW<^C2 

Denote  the  elements  of  Cn  by  an,  0  <  n  <  N  -  1,  and  characters  of  C x  by  a*, 
0<k<N-l, 

N-x  , 

ak=^2  exp{2m-^-)an . 

n~0 

Let  Djv  =  <^2,  the  dihedral  group  of  order  2N. 

Dn  =  (a,  £;  =  t2  =  1,  tat  =  a-1). 

7.1.1  Group  convolution 

For  /,  g  €  CDat, 

f=Yl  /(*)*>  9=  Y1  9ix)x »  /(*)»  $(*)  €  c, 

xGDn  xEDjv 

the  C-algebra  multiplication  is  defined  by 

f9  =  Y,  f  S  /(*)fl(*_1») )  v- 

V€Dn  \x€Dn  ) 


7.1.2  Complete  system  of  orthogonal  idempotents 

Complete  system  of  orthogonal  idempotents  of  C%  are 


(1  +  0/2,  (1-0/2. 


2  UN 


C2(ak) 


f  C2, 

{  {1}. 


k  =  0 

otherwise. 


A  complete  set  of  primitive  idempotents  of  CD  ^  are 


i{a0(l  +  0/2,  Qo(l  -  0/2,  <*n,  1  <  n  <  N  -  1}. 


^ao(l-K)  and  jftOio(l—t)  are  of  dimension  1,  while  ^an,  1  <  n  <  N  —  l 
are  of  dimension  2. 


•  N  =  2M 


k  =  0,  or  M 
otherwise. 
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Thus  complete  set  of  primitive  idempotents  of  CD2M  are 

i{a0(l  +  t)/2,  a0(l-t)/2,  aM{l  +  t)/2 , 

a^(  1  -  t)/2,an,  T<  n  <  2M  -  1,  n  ^  M}. 

27700(1  ±  £)  and  2770m (1  ±  t)  are  of  dimension  1,  while  jjan)  1  <  n  < 
N-lyn^M,  are  of  dimension  2. 

7.1.3  Basis  of  orthogonal  idempotents 

•  2  JfN 

2“o(l  +  t),an,  ia0(l  -  t),tan,  1  <  n  <  JV  -  1} 

is  a  basis  of  CD n.  The  one-dimensional  orthogonal  idempotent  basis 
elements  correspond  to  the  one-dimensional  invariant  subspaces.  The  pairs 
of  basis  elements  jf{an,tan},  1  <  n  <  N  —  1  corresponds  to  the  two- 
dimensional  invariant  subspaces. 

•  N  =  2M 

-^{^ao(l  4- 1), am,  ia0(l  -  t),tam,  ^aM (1  + 1 ), 

-  t),taM+m,  1  <  m  <  M  -  1} 

is  a  basis  of  CDw  The  one-dimensional  orthogonal  idempotent  basis 
elements  correspond  to  the  one-dimensional  invariant  subspaces.  The  pairs 
oFBasis  elements  jj{am,tam},  jj{aM+m,taM+m},  1  <  m  <  M  -  1 
corresponds  to  the  two-dimensional  invariant  subspaces. 
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Figure  1:  Profile  of  the  basis  vectors  for  C(Cz2<$  C2) 


The  matrix  of  the  new  basis  block-diagonalizes  the  group  convolution,  with 
the  computational  complexity  of  27V  log  TV.  The  diagonal  blocks  correspond  to 
the  invariant  subspaces,  and  are  of  sizes  1  or  2.  The  group  convolution  can  be 
computed  by  the  following  algorithm. 

Denote  the  matrix  of  the  new  basis  by  T(Dn).  For  f,g£  CD^,  let  h  =  fg. 

1.  Compute  /  =  r (Dw)f  and  g  =  r(D/sr)<?* 

2.  •  For  2  J(Ny  compute  h  by 

MO)  =  mm,  kN)  =  f(N)g(N), 
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h(k)  =  f(k)g(k)  +  f(N  +  k)g(2N  -  k), 

h{k  +  N)  =  f(k)g(N  +  k)  +  f(N  +  k)g(N  -  k), 

1  <  k  <  N  —  1. 

•  For  N  =  2  M,  compute  h  by 

MO)  =  mm,  h(M)  =  f(M)g(M), 

h(N)  =  f{N)g(N),  h(N  +  M)  =  f(N  +  M)g(N  +  M), 

h(k)  =  f(k)g(k)  +  f(N  +  k)g(2N  -  k), 

h(k  +  N)  =  f(k)g(N  +  k)  +  f(N  +  k)g(N  -  k), 

l<k<N-l,  ifc/M. 

3.  Compute  h  =  r(.Djv)-1h- 


7.2  The  Dihedral  Group  (Cn  x  Cn)  <)  C2 

Denote  the  elements  of  Cn  x  Cn  by  ambn,  0  <  m.  n  <  N  —  1,  and  characters 
of  Cn  x  Cn  by  ak0i,  0  <  k,  l  <  N  -  1, 


N-1N-1 

otkPi  =  7:  Y.  exp(2ni 

n— 0  m= 0 


km  4-  In 
N 


)ambn. 


(Cn-x^Cn)<$  C2  =  (a,  6,  t ;  aN  =  bN  =  t2  =  1,  a6  =  6a,  tat  =  6,  t6£  =  a). 

7.2.1  Complete  system  of  orthogonal  idempotents 

Complete  system  of  orthogonal  idempotents  of  C2  are 

(1  +  0/2,  (1-0/2. 


C2(akpl)  =  ( 


A:  =  Z,  l<kJ<N~l 
otherwise. 


A  complete  set  of  primitive  idempotents  of  C  ((C^r  x  Cn)<$  C2)  are 


^{akpk(l  + 1)/2,  akpk(  1  -  0/2,  akpu  1  <  k,  l  <  N  -  1,  k  ±  /}. 

4-  0  and  “0,  l  <  k  <  N  -  1  are  of  dimension  1,  while 

1  <  fc,  Z  <  N  —  1,  /c  ^  Z  are  of  dimension  2. 
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7.2.2  Basis  of  orthogonal  idempotents 

a0A)(l  +  t)/2,  ai 0o,  ...,aN-iPo, 
<*oPi,  aiPi{l  +  t)/2, 


<*oPn-1,  &iPn-1  •  •  • )  CtTV— l/?AT— 1(1  +  0/2. 
a0/30(l  -  f)/2,  fai/?0,  . . . ,  taN-iPo, 
ta0/3i,  ai/3i(l -t)/2,  ...,tajv_i/3i, 


tao/3N-i,  ta i(3n-i  •  •  • ,  aw-i/?AT-i(l  -  <)/2}- 

The  matrix  T  =  T  ((Cat  x  CV)<£  C2)  of  change  of  bases  is  unitary  and  has 
computational  complexity  of  4N2logN.  T  block-diagonalizes  the  group  convo¬ 
lution  into  blocks  of  sizes  1  or  2.  The  group  convolution  can  be  implemented 
by  the  following  algorithm.  For  f,gE  CD/v,  let  h  =  fg. 

1.  Compute  f  —  Tf  and  g  —  Tg. 

2.  Compute  h  by 

h(m,  m,  0)  =  /(m,  m,  0)g(m}  m,  0), 
h(m,m,  1)  =  l)g(m, m,  1),  0  <  m  <  iV  —  1, 

h(m,  n,  0)  =  /(m,  n,  0)p(m,  n,  0)  4-  /(m,  n,  l)£(n,  m,  1), 
h(m, n,  1)  sb  /(m, n,  0)£(m,  n,  1)  +  /(m,  n,  l)p(n, m,  0), 

1  <  m,  n  <  iV  —  1,  m  ^  n. 

3.  Compute  h  =  T~lh. 

7.3  (Cjv  x  Cjv)<$  C3 

G  =  (a,6,i;  aN  =  bN  =  =  1, 

ab  =  &a,£nm6n£2  =  a_nfcn”m). 

(Cjv  x  C*)*  =  {afcfl  :  0  <  fc,  l  <  N -  1}, 
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<Xk0i  =  53  exp(2nimk  +  nl)ambn. 

m,  n 

{CzY  =  {ro,  n,  r2}, 

To  =  1  +  i  +  £2 » 

n  =  1  +  wt  +  w2t2, 

r2  =  1  4-  +  wt2,  w  =  exp(27T2'/3). 

7.4  (C^<^C2)x(C^^C2) 

G  =  (a,  6, 5,  t\  aN  =  bN  =  s2  —  t2, 
ab  =  6a,  at  =  £a,  s6  =  65,  st  ~ts,  sas  =  a“2  ,tbt  =  b"1). 

Set 

<T0  =  1  +  5,  O’!  =  1  -  5  r0  =  1  +  £,  Tj  =  1  -  t. 

•  2  /at 

1- dimensional  idempotents  : 

{ao/?o^o^O)  «oj8o<To‘Ti»  <xo0q<tiTqj  aoA)0i7i} 

2- dimensional  idempotents  : 

{am/3oTo,am/30Ti,ao^n^o,ao/3n^i,  1  <  m,n  <  W  -  1}. 
4-dimensional  idempotents  : 

“  {(XmPn  :  1  <  m,  n  <  AT  -  1}. 

•  N  =  2M 

1- dimensional  idempotents: 

{ao/3o0or(b  aM/?o°oTo>  aoA)°,iro»  «M/?o^iro? 

Ofo/fo^O,  OTo,  »0/3m^170> 

ao/So^OTl,  QtAf/^OO'oTl,  O0/?0^lTl, 

<XoPm<ToTu  aM0M&OTl>  OLO0M<riTu  (*M0M<riTi}. 

2- dimensional  idempotents: 

{ot-mfl OTO?  0>  C*mA)Tl >  GLmpM'T  1 )  ^oAn^Oi  O-MfivnO 0)  1 

1  <  m  <  M  —  1,  M  +  1  <  m,  n  <  AT  —  1}. 
4-dimensional  idempotents: 

{am/3n  :  1  <  m,  n  <  M  -  1,  M  +  1  <  m,  n  <  AT  -  1}. 
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•  2  UN 

•  N  =  2M. 

Centralizers: 

(C2xC2)*,  (m,n)  =  (0,0), 

{1*}.  (m,n)  =  (0,  M),  (M,0), 

{1,  s},  {m,n)  =  {k,k), 

{1},  (m,  n)  =  (fc,  l ),  ( k ,  0),  (0,  fc), 

1  <  Jfc,  l  <  M  -  1,  M  +  1  <  k,  l  <  N  -  1. 

1- dimensional  idempotents: 

{ ql0(3q<joTo ,  QimPmVoTo,  otofloaiTo,  olmPmViTo, 

OLofio&OrU  OlmPmG  O^l,  Oiofio&lTi,  ajvf  ^McrlTl  }• 

2- dimensional  idempotents: 

{aQf3Ma0r0,  aM/30a0To,  oloPm&oTu 

{0LkPk° OR),  OtkPk&lTo}, 

{afc^_fcCToro,  ajt/^AT-AjCTiTi}, 

l<fc<  Af-1,  l  +  M<k<N-l. 

4-dimensional  idempotents: 

{otkPu  akf3o,a0f3k;  1  <  A:,  l  <  M  -  1,  1  +  M  <  fc,  /  <  N  -  1}. 

7.6  (C^v  x  Cjv)<^  ^4 

G  =  (a,b,t-,  aN  =  bN  =  t4, 

ab  =  ba,tambnt3  =  anb~m). 

8  Results  of  numerical  experiments 

The  following  pictures  illustrate  the  power  of  our  methods  to  discriminate  be¬ 
tween  object  and  noise  in  image  data.  Generally,  we  expect  that  image  infor¬ 
mation  will  distribute  over  direct  sum  factors  with  noise  more  prominent  in  low 
dimensional  factors  and  object  more  concentrated  in  high  dimensional  factors. 
In  all  but  one  of  the  illustrations  the  object  can  be  recovered  from  the  highest 
dimensional  factor  with  significant  noise  reduction  as  measured  by  the  signal  to 
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noise  ratio  (SNR).  The  one  in  which  this  is  not  the  case,  the  second  highest  di¬ 
mensional  factor  becomes  most  significant.  The  anomaly  reflects  the  directional 
bias  of  the  imaging  model  and  emphasizes  the  importance  of  the  imaging  model 
as  a  tool  for  highlighting  certain  image  characteristics. 

The  two-dimensional  images  consists  of  delta  functions,  lines  of  delta  func¬ 
tions  and  crossing  lines  of  delta  functions  embedded  in  noise.  Delta  functions 
have  pixel  support  and  should  be  distinguished  from  regional  boundaries.  The 
noise  is  at  30%  variance  to  the  object.  The  imaging  model  indexes  each  coor¬ 
dinate  direction  by  the  dihedral  group  Cg  <$  C2  producing  one,  two  and  four- 
dimensional  irreducible  factors. 
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G  =  (C8^C2)x(Cs4C2) 

Figure  2:  Images  in  the  following  are  1)  Single  pixel  of  value  1,  embedded  in 
noise  of  values  0-0.25;  2)  1-dimensional  irreducible  factors;  3)  2-dimensional 
irreducible  factors  and  4)  4-dimensional  irreducible  factors. 
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C  =  (C8^C2)  X  {C^C2) 

Figure  3:  Images  in  the  following  are  1)  Two  pixels  of  value  1  and  1.5,  embedded 
in  noise  of  values  0-0.5;  2)  1-dimensional  irreducible  factors;  3)  2-dimensional 
irreducible  factors  and  4)  4-dimensional  irreducible  factors. 
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G  =  (C84C2)x(C84C2 ) 

Figure  4:  Images  in  the  following  are  1)  collection  of  pixels  of  value  1  embedded 
in  noise  of  values  0-0.3;  2)  1-dimensional  irreducible  factors;  3)  2-dimensional 
irreducible  factors  and  4)  4-dimensional  irreducible  factors. 


G  =  (C8^C2)x(Ca^C2) 

Figure  5:  Images  in  the  following  are  1)  collection  of  pixels  of  value  1  embedded 
in  noise  of  values  0-0.3;  2)  1-dimensional  irreducible  factors;  3)  2-dimensional 
irreducible  factors  and  4)  4-dimensional  irreducible  factors. 


G=(C8<^C2)x(C8^C2) 

Figure  6:  Images  in  the  following  are  1)  collection  of  pixels  of  value  1  embedded 
in  noise  of  values  0-0.3;  2)  1-dimensional  irreducible  factors;  3)  2-dimensional 
irreducible  factors  and  4)  4-dimensional  irreducible  factors. 
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The  following  pictures  illustrate  the  power  of  our  methods  as  a  data  reduc¬ 
tion  tool  for  textures.  Generally,  textural  information  concentrates  in  low  di¬ 
mensional  factors  with  higher  dimensional  factors  representing  textural  changes. 
The  bias  of  the  imaging  model  affects  which  factors  concentrate  information  but 
not  the  number  of  coefficients  required  to  encode  the  information.  Generally 
by  slightly  modifying  the  imaging  model  or  by  repositioning  the  texture,  most 
information  upto  some  resolution  can  be  placed  in  the  lowest  dimensional  factor. 
Since  the  number  of  coefficients  required  to  encode  this  factor  is  a  small  fraction 
of  the  number  of  coefficients  describing  the  texture,  significant  data  reduction 
is  possible.  In  Figure  12,  we  display  the  16  coefficients  required  to  describe  a 
32  x  32  coefficient  texture. 
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G  =  {C8$  C2)  X  (C8  4  c2) 

Figure  7:  Images  in  the  quadrants  are  respectively  1)  Simulated  image  of  values 
0,  1,  2  and  3;  2)  1-dimensional  irreducible  factors;  3)  2-dimensional  irreducible 
factors  and  4)  4-dimensional  irreducible  factors. 


68 


G  =  (C84C2)x(Cs4C2) 

Figure  8:  Images  in  the  quadrants  are  respectively  1)  Simulated  image  of  values 
0-6;  2)  1-dimensional  irreducible  factors;  3)  2-dimensional  irreducible  factors 
and  4)  4-dimensional  irreducible  factors. 
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G  =  (C8<*C2)x(C8^C2) 

Figure  9:  Images  in  the  following  are  1)  Simulated  image  of  values  0-6  de¬ 
scribing  texture  with  noise  values  of  0-.7  ;  2)  1-dimensional  irreducible  fac¬ 
tors;  3)  2-dimensional  irreducible  factors  4)  4-dimensional  irreducible  factors; 
5)  2-dimensional  irreducible  factors  uniformly  multiplied  by  3;  6)  4-dimensional 
irreducible  factors  multiplied  uniformly  multiplied  by  10.  , 
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G=(CHC2)x(CHC2) 

Figure  10:  Images  in  the  following  are  1)  Simulated  image  of  values  0,  1,  2 
describing  texture  with  noise  values  of  0-0.25  ;  2)  1-dimensional  irreducible 
factors;  3)  2-dimensional  irreducible  factors  and  4)  4-dimensional  irreducible 
factors. 
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G  =  (C8^C2)x(C8<$C2) 

Figure  11:  Images  in  the  following  are  1)  Simulated  image  of  values  0,  1,  2 
describing  texture  with  noise  values  of  0-0.25  ;  2)  1-dimensional  irreducible 
factors;  3)  2-dimensional  irreducible  factors  and  4)  4-dimensional  irreducible 
factors. 
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G  =  (C8<*C2)x(C8<*C2) 

Figure  12:  Images  in  the  following  are  1)  Simulated  image  of  value  10  describing 
a  square  with  noise  values  of  0-0.25  ;  2)  1-dimensional  irreducible  factors;  3) 
2-dimensional  irreducible  factors  and  4)  4-dimensional  irreducible  factors. 
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G  =  (C8^C2)x(Cs4C2) 

Figure  13:  Images  in  the  following  are  1)  Simulated  image  of  value  10  describing 
a  square  with  noise  values  of  0-0.25  ;  2)  1-dimensional  irreducible  factors;  3) 
2-dimensional  irreducible  factors  and  4)  4-dimensional  irreducible  factors. 
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G^(C8^C2)x(Cs4C2) 

Figure  14:  Images  in  the  following  are  1)  Simulated  image  of  value  10  describing 
a  square  with  noise  values  of  0-0.25  ;  2)  1-dimensional  irreducible  factors;  3) 
2-dimensional  irreducible  factors  and  4)  4-dimensional  irreducible  factors. 
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G  =  (C8^C2)  X  (Cs4C2) 

Figure  15:  Images  in  the  following  are  1)  Simulated  image  of  value  10  describing 
a  square  with  noise  values  of  0-0.25  ;  2)  1-dimensional  irreducible  factors;  3) 
2-dimensional  irreducible  factors  and  4)  4-dimensional  irreducible  factors. 


The  following  two  figures  illustrate  the  power  of  filtering  in  the  coefficients 
of  image  expansions  relative  to  the  basis  of  translated  idempotents  to  data 
compression  and  noise  reduction. 


G  =  (Cio  4  C2)  x  (C10  <£  C2) 

Figure  16:  Images  in  the  following  are  1)  Simulated  image  of  value  10  describ¬ 
ing  a  square  with  noise  values  of  0-1  ;  2)  Non-zeros  coefficients  describing  the 
image;  3)  Reconstruction  using  the  100  largest  coefficients  in  the  two  larger  di¬ 
mensional  irreducible  factors;  4)  Location  of  the  non-zeros  coefficients  used  in 
the  reconstruction. 
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G=(C10^C2)x(C10^C2) 

Figure  17:  Images  in  the  following  are  1)  Simulated  image  of  values  of  5  and 
10  describing  a  square  with  noise  values  of  0-1  ;  2)  Non-zeros  coefficients  de¬ 
scribing  the  image;  3)  Reconstruction  using  the  8  largest  coefficients  in  the  one¬ 
dimensional  irreducible  factors;  4)  Location  of  the  non-zeros  coefficients  used  in 
the  reconstruction. 
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